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Preface 



Few-body problem plays a very important role in the physics of ultra-cold gases. 
It describes the structure of molecules and their formation in three-body collisions 
(three-body recombination), atom-molecule and molecule- molecule collisional proper- 
ties, structure of trimers and larger clusters, numerous phenomena which go under the 
name of Efimov physics, and many other problems. 

In ultracold gases we benefit from a remarkable separation of scales. The atomic 
do Broglie wavelengths are much larger than the range, R e , of the van dcr Waals 
interatomic potential. The ultracold regime is defined by the inequality kR e <C 1, where 
k is the typical atomic momentum. In this regime very few terms in the effective range 
expansion for the scattering amplitude suffice, and most of the time the interaction 
is characterized by a single parameter - the s-wave scattering length a. From the 
two-body viewpoint all short-range potentials are equivalent as long as they have the 
same scattering length, and, therefore, one can use an idealized zero-range potential 
(pseudopotcntial) with the same a. The zero-range model well describes the weakly 
interacting BEC (Huang 1963) as well as the whole range of the BCS-BEC crossover 
in fermionic mixtures (Giorgini et al. 2008). 

There is a class of few-body problems where both long and short lengthscales 
are important. For example, the knowledge of a is not sufficient for calculating the 
spectrum of Efimov trimers or the rates of recombination and relaxation to deeply 
bound molecular states. Such problems can be solved in the universal limit, \a\ 3> R e , 
by introducing the so-called three-body parameter, which absorbs all the short-range 
three-body physics in the same manner as the scattering length absorbs the short- 
range two-body physics. Universality (Braaten and Hammer 2006, 2007) in this context 
reflects the amazing fact that different systems with the same scattering length and 
three-body parameter exhibit the same physics. The possibility to modify a in atomic 
gases by using Feshbach resonances makes ultracold gases an ideal playground to check 
the universal theory. The Efimov effect predicted 40 years ago (Efimov 1970) has been 
first observed in a cold gas of 133 Cs in Innsbruck (Kraemer et al. 2006, see also Ferlaino 
et al. 2011 for review) and subsequently in other alkali atoms and mixtures (Ottenstein 
et al. 2008; Huckans et al. 2009; Zaccanti et al. 2009; Barontini et al. 2009; Gross et 
al. 2009, 2010; Pollack et al. 2009; Nakajima et al. 2010, 2011; Lompe et al. 2010). 

The few-body analysis, apart from being interesting on its own, can be used in 
many-body problems to integrate out few-body degrees of freedom, thus making the 
many-body problem more tractable. For example, a two-component fermionic mixture 
on the BEC side of the BCS-BEC crossover is actually a Bose gas of molecules (or an 
atom- molecule mixture in the density imbalanced case). If the molecules are sufficiently 
small (a is much smaller than the interparticle distance), we can forget about their 
composite nature and treat them as elementary objects as long as the atom-molecule 
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and molecule-molecule scattering parameters are known, the latter being given by the 
solution of the corresponding three- and four-atom problems. Here we profit from the 
fact that the few- and many-body processes separate on the length and energy scales. 
Then the equation of state and the stability of the system with respect to collapse 
or phase separation are much easier to analyze. We can also mention lattice problems 
where in order to describe the tight binding limit one has to solve the few-body problem 
on a single site and use the result as an input for the Hubbard model. 

Another path from "few" to "many" is the high-temperature virial expansion, 
which consists of solving few-body problems with more and more particles (Huang 
1963; Liu et al. 2010a). The results are then used for consecutive approximations 
of thermodynamic quantities. The procedure rapidly converges, which emphasizes the 
importance of the few-body aspect of many-body problems and leads to the conjecture 
that even at low temperatures and in the strongly interacting regime, where we lack a 
small parameter, the natural sorting criterion for different sets of Fcynman diagrams 
is the number of particles (particle-hole pairs) involved. 
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1 



The two-body problem and 
resonance width 



In this chapter we discuss the effective range expansion for the two-body scattering 
amplitude and put a special emphasis on the role of the resonance width. Near a scat- 
tering resonance the scattering length can be modified and, in particular, can take 
anomalously large values (i.e. \a\ 3> R e )- Another very important parameter charac- 
terizing the resonance is its width, which is determined by the strength of the coupling 
between the closed and open channels. The narrower the resonance, the stronger the 
collision energy dependence of the scattering amplitude. In the effective range expan- 
sion this leads to an increased value of the effective range. 

In order to familiarize ourselves with the notions of scattering length, resonance 
width, and effective range we will consider a simple model of a single-channel potential 
with a barrier. It mimics the Feshbach resonance picture and clearly illustrates the role 
of the barrier. Then we will derive a general formula for the effective range and discuss 
the relation between the effective range and the probability of finding two atoms in the 
closed channel. This quantity plays a very important role in many problems involving 
Feshbach resonances with finite width. As a particular example, we will discuss the 
structure of a Feshbach molecule. 

1.1 Potential well with delta-function barrier 

Consider two atoms with reduced mass fj, interacting via an isotropic interaction po- 
tential V(r). The radial Schrddingcr equation for the relative motion with zero orbital 
angular momentum reads 



1 B 2 PI 



(i.i) 



where E = k 2 /2fx is the collision energy and we use x(r) = rip(r). 
Let us take V(r) in the form 



V(r)=g6(r-R e )-^ 




(1.2) 



The potential (|1.2p is plotted in Fig. 11.11 It allows one to vary the resonance width at 
will by varying the barrier strength g. The solutions of Eq. (|1.1[) in regions I and II 
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V(r)k 




Fig. 1.1 The square well potential with a delta- function barrier (|1.2[) and the wavefunction 
in the inner, xi, an d the outer, xn, regions at the collision energy fc 2 /2/i (solid line). The 
spike of the wavefunction at r = R e is governed by Eq. (|1.3p . The dashed line shows the 
wavefunction xn extrapolated to the inner region I. 



are, respectively, xi ^ sin(-y/ Kq + k 2 r) and xn oc sin(fcr + 8). The matching condition 
at r = R e reads 

X ^(Re)/Xll(Re) - x'l(Re)/Xl(Re) = 9, (1.3) 

which gives the phase shift S: 



cot (5 = 



\l + k 2 + tan(v/K^ + k 2 R e )[g + k tan(fci? e )] 
k — tan(kR e )[g + \/ + k 2 tan(y / Kp + k 2 R e )} 



(1.4) 



The s-wave scattering amplitude is given in terms of the phase shift as (Landau and 
Lifshitz 1987) 

f(k) = l/(kcot6(k)-ik). (1.5) 

The idea of the zero-range approximation is to extrapolate the solution xn into 
region I (dashed line in Fig. II. ip and use its logarithmic derivative at the origin to 
construct the zero-range boundary condition which then replaces the potential V: 



X'(0)/X(0) = k cot 6(k). 



(1.6) 



By construction, the s-wave scattering amplitudes obtained by solving Eq. (jl.ll) on 
the one hand and the free-motion Schrodingcr equation with the boundary condition 
(jl.6| on the other hand are identical. In general, substituting the potential by the 
boundary condition (|1.6j) docs not make the scattering problem easier since we do 
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have to calculate the phase shift. However, the zero-range approximation becomes 
valuable for small momenta, when one needs only a few parameters to describe the 
scattering. Namely, one writes down the effective range expansion 

kcotS(k) = -l/a+ (r /2)fc 2 + (1.7) 

where a is the scattering length, ro is the effective range, and the terms denoted by ... 
contain higher powers of k. 

Typically (i.e. for most of interatomic interactions provided by nature), the length 
parameters a and ro are of the order of the physical range of the potential R e . Then, 
in the ultracold limit the first term on the right hand side of Eq. (|1.7p suffices and one 
arrives at the so-called Bethe-Peierls boundary condition for the wavefunction 

x'(o)/x(o) = N>M)7^Mk=o = -i/o, (1.8) 

or, equivalcntly, 

ip(r) oc l/r- 1/a, r ->■ 0. (1.9) 

Equation (|1.9p sets the relation between the coefficient in front of the singular term 
and the regular term of the wavefunction. Besides, it points to the physical meaning 
of the scattering length - for positive a the node of the zero-energy wavefunction is 
found at r = qj. In this sense, at sufficiently low energies any potential characterized 
by a > is equivalent to a hard-core one of radius a. 

In the case of the potential ()1.2[) the scattering length is given by 

a = R e - l/[g+ K cot(n R e )]. (1.10) 

In Fig. 11.21 we show the dependence a(«o) for g = (left panel) and gR e = 10 (right 
panel). We see that by changing the depth of the well we can tune the scattering length 
to a resonance (a = oo). This happens every time a new bound state crosses the zero 
energy threshold. 

In the limit g = oo the regions I and II completely decouple, the former we can 
call closed channel, the latter - open channel. Obviously, in this case stationary states 
in the closed channel live independently and the potential (jl.2[) is equivalent to a 
hard-core one characterized by a = R e . For finite but large g the stationary states 
at positive energies become quasi-stationary and when the collision energy coincides 
with the position of one of them we encounter a narrow scattering resonance. This 
type of resonant scattering is called the Breit-Wigner scattering or the scattering on 
a quasi-stationary level. 

Expanding Eq. (|1.4p at small momenta one obtains the effective range ro for our 
model potential ()1.2|) . In fact, it makes sense to discuss ro only for very large values 
of a when 1/ak is comparable to r$k [see Eq. (|1.7j) ]. In such a narrow vicinity of the 
resonance ro can be considered equal to its value at a = oo: 

r (a = oo) = R e - g(l + gR e )/n 2 ► -(g / k ) 2 R e - (1.11) 

g->oo 

^^Here we mean the wavefunction obtained by extrapolating its large-distance asymptote to the 
region r ~ R e . 
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Fig. 1.2 The scattering length versus no for g = (left panel) and gR e = 10 (right panel). 



We see that in the limit of a narrow resonance (g — > oo) the effective range ro is 
negative and can be much larger than the physical range of the potential. Below we 
will see that this statement holds generally for narrow resonances. 
Let us introduce a positive length parameter 



R* = -r /2 > 0. 



(1.12) 



We call a resonance narrow if R* ^> R e . Keeping only the first two terms in the 
effective range expansion (|1.7[) we write the scattering amplitude (|1.5[) in the form 



/(*) = - 



1 



1/a + R*k 2 + ik' 



(1.13) 



Equation (|1.13l) has the well-known low-energy Breit-Wigner resonance shape, which 
becomes apparent if we rewrite it equivalently in energy units: 



HE) = 



E-E r 



ijyE 



(1.14) 



where the position of the quasi-stationary state E res and the tunneling amplitude 7 
are related to the lengths scales a and R* as a = —hj/ y/2jlE res and R* = h/^/2jl-f. 
The parameter r yy/~E (multiplied by 2) is nothing else than the decay rate of the quasi- 
stationary state, and therefore, both 7 and R* are necessarily positive. 

Similar to resonances in the model potential with barrier (|1.2[) magnetic Feshbach 
resonances occur when the collision energy of two atoms is close to the energy of a 
quasi-discrete molecular state in another hypcrhnc domain, in this case closed channel. 
The tuning of the scattering length is achieved by shifting the open and closed channels 
with respect to each other in an external magnetic field (hypcrhnc states corresponding 
to the open and closed channels have different magnetic moments). For magnetic 
resonances 

R * = 9 ( L15 ) 

where a bg is the background scattering length, /i rol is the difference in the magnetic 
moments of the closed and open channels, and AB is the magnetic width of the 
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Fcshbach resonance. Usually, there are many Feshbach resonances for a given atomic 
gas or mixture in the realistic magnetic field range and, most of the time, one has 
a possibility to work with a couple of "nice" wide resonances. However, there are 
systems where all available resonances are narrow. For example, parameters of the 
widest resonance in 87 Rb (at 1007.4 G) arc A B w 0.17 G, R e « 4.4nm, R* w 32nm 
(Volz et al. 2003). Similarly, all 6 Li- 40 K interspecies resonances discussed so far are 
characterized by R* > lOOnm (Wille et al. 2008; Tiecke et al. 2009), whereas the 
van der Waals range for this alkali pair is R e w 2.2nm. The review of Chin et al. 
(2010) broadly covers the phenomenon of Feshbach resonances and related physics in 
ultracold gases and, in particular, presents parameters for many resonances. 

As we have already mentioned, the situation where a is much larger than the phys- 
ical range of the potential is rare and it usually requires a finctuning of the potential. 
Large R* arises for narrow resonances when the scattering amplitude is character- 
ized by a strong dependence on the collision energy. One can imagine an interatomic 
potential for which the higher order terms in the effective range expansion (|1.7jl arc 
also anomalously large. For example, we can introduce one or several additional closed 
channels with quasi-stationary states very close to the open channel threshold result- 
ing in a rather exotic scattering amplitude. However, in all practically relevant cases 
these terms can be neglected as the typical distance between neighboring molecular 
states is rather large - of order \j [xR\. In particular, for the model potential ()1.2[) the 
next term in the expansion ()1.7[) at a = oo and in the limit of large g is proportional 
to (5 l //to) 4 i?gfc 4 cx R* 2 R e k 4 . This means that the potential (|1.2j) has a well defined 
zero-range limit: tending R e to zero one can accordingly modify kq and g in such a 
way that a and R* remain unchanged whereas all other terms in the effective range ex- 
pansion vanish. In this limit Eqs. (|1.13|) and (|1.14[) become exact and the Bethc-Peierls 
boundary condition reads 

X'(0)/X(0) = (r^(r))7nKr)| r=0 = -1/a - 2^R*E, (1.16) 

or, in analogy with (|1.9I) . 

V>(r) cx 1/r- 1/a- 2fj,R*E, r -> 0, (1.17) 

where E is the collision energy. 

Clearly, a is not always enough to characterize the interaction strength. For ex- 
ample, a can be infinite, but typical energies in the problem can be so high that the 
interactions are off-resonant. Thus, as a more appropriate quantity it is reasonable to 
introduce the so-called effective energy-dependent scattering length a(E) defined by 
the equation 

l/a(E) = 1/a + 2[iR*E. (1.18) 

Comparison of a(E) with other lengthscalcs in the problem gives us an idea of how 
strong the interaction is at a given (or typical for this problem) collision energy. 

1.2 Effective range and population of closed channel 

Let us now demonstrate the relation between the effective range and the probability 
to find two atoms within the physical range of the potential. We will then use this 
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relation to understand the structure of a weakly bound molecular state near a narrow 
resonance. 

Consider the Schrodingcr equation with an arbitrary short-range potential 

V(r). We require that V(r) can be neglected at r > R e and that the wavefunction \ 
vanishes at the origin. Assume that xo and Xk are two real solutions of Eq. (|l.ip at 
zero and finite energy, respectively: 

[-d 2 /dr 2 + 2 f iV(r)} X o(r) = 0, (1.19) 
[-d 2 /dr 2 + 2fiV(r) - k 2 ] Xk (r) = 0. (1.20) 
We multiply Eq. (|1.19j) by Xk and Eq. (|1.20[) by xo- Subtracting the results we obtain 

(-Xox'k+XkXo)' = k 2 X oXk, (1.21) 

Let us now introduce the functions xo and Xk which satisfy, respectively, Eqs. (|1.19l) 
and (|1.20|) with V(r) = and which equal xo and Xk at distances r > R e . These 
functions satisfy 

(-Xox'k + XkXoY = k 2 X0Xk- (1-22) 

We subtract Eq. (|1 .22[) from Eq. (|1.21[) and integrating the result from to R e we 
obtain 

f R e 

-Xo(0)x' fc (0)+x fc (0)x'o(0)-fc 2 / [xo(r)Xk(r) - Xo(r)Xk(r)]dr. (1.23) 



In the limit k — > we can use the boundary condition ([1.6)1 for xo and Xk and the 
effective range expansion ([1.7)1 up to the effective range term. For the effective range 
we obtain 

R* = -ro/2 = Xo 2 (0) ^ 6 Xl(r)dr - xl(r)dr\ . (1.24) 

The first integral on the right hand side of Eq. (|1.24p gives the probability of finding 
the two atoms in the region r < R e : 



r<R c 



Re 

2 



in / xt{r)dr. (1.25) 



Near a narrow resonance P r <R e is dominated by the population of the quasi-stationary 
state in the closed channel. The second integral in Eq. (jl.24[) equals 



Re 

Xl(r)dr = Xo(0)Re[l - R e /a + (i? e /a) 2 /3], (1.26) 

where we have used the equality xo( r ) = Xo(0)(l — r/a). In the case R* 3> R e the 
second integral can be neglected and we get 

P r <R c ~^xl(0)R*. (1.27) 

Equation (jl.27[) allows one to find the occupation of the closed channel (sometimes 
called bare molecular state) without solving coupled channel equations. One just needs 
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to know the effective range parameter R* , which can be found from Eq. (|1.15[) , and the 
short distance asymptote of the open channel wavefunction, ip(r) « xo(0)(l/r — 1/a). 
Note, that Xo(0) is nothing else than the coefficient in front of the singular part of the 
wavefunction of the two-particle relative motion just outside of the potential support . 

Equation (|1.27j) becomes exact in the zero-range limit and, together with Eqs. (|1.13p . 
([1.140 . and (|1.16p it generalizes the usual zero-range theory to the case of narrow 
resonances. Let us now demonstrate how this modified zero-range theory works by 
calculating parameters of a weakly bound molecular state near a narrow resonance. 

1.3 Weakly bound molecule near a narrow resonance 

The wavefunction and the energy of a weakly bound molecular state can be found in 
the zero-range approximation, the applicability of which in this case requires that the 
size of the bound state be much larger than the physical range of the potential R e . 
The solution of the free-motion Schrodinger equation 

[~d 2 /dr 2 + K 2 } X (r) = 0, (1.28) 

decaying at large distances reads 

X(r) = a norm exp(-Kr), (1-29) 

where a norm is a real normalization coefficient and k is determined from the Bcthe- 
Pcicrls boundary condition ([1.160 : 

x'(0)/x(0) = - K = -l/a + R* K 2 . (1.30) 

This quadratic equation can also be obtained by looking for the poles of the scattering 
amplitude 1)1.140 . One of its roots gives the energy of the bound state 

e = -^=-±(£±M^l)\ (L31) 
2fi 2/i ^ 2R* I V ' 

The normalization constant a is calculated in the following way. The integral of 
the square of x gives the probability of finding the atoms in the "open" channel 
-Popen = Pr>R e (i n this case R e = 0): 

/>oo 

Popen = 4tt / x 2 (r)dr = 27rc£ omi //s. (1-32) 



At the same time the closed channel population is given by Eq. (|1.27|) : 

Closed = 47r X 2 (0)i?* = A-Kal OTm R*. (1.33) 
Using Eq. (11.311) and the fact that P op cn + ^closed = 1 we get 



Popen = l/y/l+4R*/a, (1.34) 

and 



1 - l/Wl + AR*/a 



(1.35) 



2tt(1 + 2R*k) V ^R* 
Depending on the value of the ratio R* /a we distinguish two regimes: the case 
R* /a <C 1 we call the regime of small detuning, and the opposite limit - the regime 
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of intermediate detuning. To be more precise, in the latter case we require R e <C 
a <C R* since we would like to stay within the region of validity of the zero-range 
approximation. In the regime of small detuning the molecular energy equals £o ~ 
— l/2fia 2 and the probability of finding the atoms in the closed channel is small. The 
physics in this case is similar to the one near a wide resonance. In the regime of 
intermediate detuning £o ~ — 1/2 fiR* a Rj E res and P c iosed ~ 1, which means that the 
molecular state in this case is dominated by the closed channel or bare molecular state. 



2 



Basics of the three-body problem 
with short-range interactions 



In this chapter we apply the zero-range approach to the problem of three atoms. We 
consider the phenomenon of the three-body collapse, or Thomas collapse, and define 
the Efimov effect. Then we introduce a model system of one light and two heavy 
atoms where one can use the Born-Oppcnhcimcr approximation. The Efimov physics 
in this model is extremely transparent and one can understand the role of the two-body 
interaction parameters (scattering length, resonance width), masses of the atoms, and 
their quantum statistics. We also discuss the question of inelastic three-body losses. 

2.1 Thomas collapse 

Consider three identical bosons interacting with each other via the square well poten- 
tial (|1.2p with .9 = 0. Let us gradually decrease R e and at the same time increase the 
depth of the potential in such a way that the scattering length a is unchanged. In this 
way we arrive at a zero-range potential equivalent to the initial one in the ultracold 
limit. What happens with the spectrum of the three-body problem as we decrease 
R e l Below we show that its ground state energy is unbound from below in the limit 
R e — > 0. This phenomenon was discovered by Thomas (1935) who theoretically ad- 
dressed the mass defect of tritium and arrived at a number of important conclusions 
on the character of the proton-neutron and neutron-neutron interaction. 

The Hamiltonian for three bosons in the center-of-mass reference frame reads (we 
set h = m = 1) 

H = -V* - V 2 + V(y) + £ V(\V3x ± y|/2), (2.1) 

± 

where y is the distance between two of them and -\/3x/2 is the distance from their 
center of mass to the third boson. Consider the variational wavefunction 

*var(x,y) = (w/2^) 3 / 2 cxph W (a; 2 + y 2 )/4], (2.2) 

which is the normalized ground state wavefunction of a six-dimensional harmonic 
oscillator of mass 1/2 and frequency lo - our variational parameter. The ground state 
energy of the Hamiltonian (|2.1[) satisfies 

E < (tfvarlflltfvar) = 3w/2 + 3(w/2^) 3 / 2 ^ V(y) cxp(-ujy 2 /2)d 3 y 

= 3w/2 - 3fi^[erVw-R 2 /2 - v / 2wi? 2 /7rexp(-wi? 2 /2)], (2.3) 
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where we use Eq. (|1.2[) and the fact that the three interaction terms in (|2.1j) contribute 
equally to the matrix element (^ var \H l^var)- 

When we reach the condition R e <C \a\ the depth of the potential becomes practi- 
cally independent of the exact value of a and is determined by the equation cot(Koi?e) = 
[see Eq. (|1.10[) ]. The minimum depth of the potential is then given by kq = n/2R e . 
Substituting this value to Eq. (|2.3j) we find that its right hand side is minimized for 
w « 2.1/Rl and we get 

E Q < -0.167/Rl, (2.4) 

which confirms the statement that the ground state energy is unbound from below as 
R e -> 0. 

In the case of three bosons this result is expected: The inequality \a\ ^> R e means 
that there is a bound or quasi-bound two-body state. In 3D, this requires that the 
potential is of a certain finite attractive strength. If we a priori assume that two of the 
atoms are very close to each other, the third atom will experience twice the attractive 
potential and will then be tightly bound, which verifies the initial assumption. In- 
terestingly, Thomas showed that the statement holds for the proton-neutron-neutron 
(tritium) system in which the neutron-neutron interaction is neglected. In this case the 
above hand waving argument does not work because the interaction potential in be- 
tween a neutron and a deutron is the same as the neutron-proton interaction. However, 
the Thomas effect can be explained by the larger effective mass of the neutron-deutron 
system compared to the neutron-proton one. The strict variational analysis is rather 
bulky in this case (Thomas 1935). 

The system of three 4 He atoms is one of the most spectacular examples in which 
this physics actually takes place. The 4 He- 4 He potential is just attractive enough to 
accommodate a single bound state close to the continuum, the scattering length is 
positive and is an order of magnitude larger than the range of the potential. However, 
the 4 He3 trimer is much smaller in size and significantly stronger bound than the 
dimer. In fact, even if the interatomic potential were a bit weaker, so that the dimer 
state is just pushed into the continuum, the trimer state would still be there. This 
phenomenon is called the Borromean binding - none of the two-body subsystems has 
a bound state, but there exists a trimer state. An example of the Borromean binding 
is the 6 He nuclei which is stable in spite of the fact that neither a-n nor n-n system 
is bound. 

We now know that the ground state of such systems is a trimer state with the 
binding energy ~ \jR\. Are there other trimer states and what does their spectrum 
look like? 

2.2 Efimov effect and discrete scaling invariance 

Consider again three identical bosons interacting with each other via a potential of 
the range R e . This time let us keep R e constant, but increase the value of a, no 
matter positive or negative, by slightly tuning the depth of the potential. Efimov 
showed (Efimov 1970) that as a — > oo the number of trimer states is proportional 
to (so/7r) log(|a|/i? e ) — s- oo, where sq is a dimcnsionlcss parameter. A new trimer 
state emerges from the continuum each time a is increased by a factor of exp(7r/so). 
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When a = oo we have infinite number of trimer states with binding energies \e n \ oc 
cxp(— 2im/ sq). 

The effect is quite general - it can take place in a system of three particles with 
different masses and statistics as long as at least two of the three interactions are 
resonant. The parameter so depends on the masses and statistics of the atoms, for 
identical bosons it equals s ~ 1-00624 so that the scaling parameter exp(7r/s ) ~ 22.7. 

In fact, the Efimov effect goes along with the phenomenon of discrete scaling in- 
variance: a three-body observable of the dimension of length to the power a close to 
the resonance (a — > oo) behaves as a a f(a/R e ), where / is a dimensionless function 
periodic on the logarithmic scale, i.e. f(x) = f(xexp[ir/so\). 

We are used to the idea of continuous scaling invariance: once we multiply all 
lengthscales in the Hamiltonian by a factor of A, all eigenenergies get multiplied by 
1/A 2 and all eigenfunctions trivially rescale to satisfy the normalization condition. It 
is important to keep in mind that in the example discussed above we do change a, but 
we do not rescale R e . Therefore, the discrete scaling invariance does not contradict 
our understanding of the continuous scaling. As we show below, when R e —¥ the 
short-range physics does not completely drop out of the (three-body) problem. 

2.3 Born-Oppenheimer approximation 

Let us consider a light atom of mass m interacting resonantly via a short-range po- 
tential with two heavy atoms of mass M. For M/m ^ 1 one can use the Born- 
Oppenheimer approximation to illustrate the essential physics leading to the Efimov 
effect (Fonseca et ai. 1979). The Born-Oppenheimer approach takes advantage of the 
large mass ratio by assuming that the state of the light atom adiabatically adjusts 
itself to the distance R between the heavy atoms. The three-body problem thus splits 
into two "one-body" pieces. First, we solve the problem of the light atom in the field of 
two fixed scatterers at positions R/2 and — R/2 (a version of the double- well potential 
problem) . The binding energy of the light atom then serves as the interaction potential 
between the heavy ones. The second step is then to solve the Schrodinger equation for 
the relative motion of the heavy atoms. 

If there were only one heavy atom at position R/2, the wavefunction of the light 
atom could be written as 

-re|r-R/2| 

(2 - 5) 

where k satisfies Eq. (|1.30[) . In the case of two heavy atoms we look for the solution 
in the form 

p -«(ft)|r-R/2| -re(fl)|r+R/2| 

« Cl |r-R/2| + ° 2 |r + R/2| ' ™ 
The wavefunction (|2.6[) satisfies the Schrodinger equation 



■^ + £mr±R/2| 



V'rW = e(R)^ R (r) (2.7) 



with the energy 
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e(R) = -n 2 (R)/2m. (2.8) 

The parameters n(R), Ci, and C2 are obtained by applying the Bethe-Peierls boundary 
condition (|1.16j) to ipn(r) at vanishing r ± R/2 and by using the fact that the collisional 
energy equals e(R). In particular, for the boundary r — R/2 — > the resulting equation 
reads 

[k(R) - l/a+ R*K 2 (R)]d - {exp[-K(R)R]/R]}C 2 = 0, (2.9) 

and for the other boundary we get the same equation with interchanged C\ and C2 . 

These equations have two solutions: {Ci, 6*2}+ = {1, 1} and {Ci, C2}- = {1, — 1}. 
Substituting them into Eq. (|2.6[) we obtain the wavefunctions -0R i+ (r) and ?/>R,_(r), 
which are nothing else than the symmetric and antisymmetric bound states of our 
double- well potential problem. The equation for n± (R) reads 

n ± (R) t cxp[-K ± (R)R] /R = 1/a- R*k±(R). (2.10) 

Solving Eq. (|2 . 10[) and using Eq. (|2 .8[) we get the induced interaction potentials e± (R) 
for the heavy atoms. In Fig. 12.11 we plot these potentials for different values of R* /a. 
We observe that e+(R) is purely attractive and e~{R) is repulsive. In the limit of 
large R both potentials tend to the molecular binding energy (note that the Born- 
Oppenheimcr value of eo differs from Eq. (|1.31[) by the factor fj,/m). At distances 
R < a the symmetric state (+) is the only bound state since the antisymmetric one 
(-) goes to the continuum at R = a. 

1 




af -1 

5 



-3 
-4 

12 3 

R /a 

Fig. 2.1 Born-Oppenheimer potentials e_ and e+ in units of the molecular binding energy 
I eo I for different detunings: R* /a — (solid), R* /a = 1 (dashed), and R* /a — 10 (dotted). 

The second step of the Born-Oppenheimer method consists of solving the Schrodingcr 
equation for the heavy atoms by using e± (R) as the potential energy surface. Let us 
postpone the question on how this induced interaction depends on R* /a and first con- 
sider the case R* = 0. Then, in the limit R^> a the parameter k±(R) can be expanded 
in a power series in the small parameter exp(— R/a). The first three terms read 
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K±(R) « a" 1 ± JT 1 exp(-R/a) - R' 1 cxp(-2i?/a), (2.11) 
which gives the long-distance asymptotes: 

f <lt\~ 1 _ exp(-i?/q) exp(-2g/a) 

£±(i?) ~ ~2W mai? + m aR (1 ~ a/2i?) ' (2J2) 

This means that at large distances the exchange of the light atom leads to a Yukawa- 
type force between the heavy atoms. In the limit fi«a only the symmetric state is 
bound, the right hand side of Eq. (|2.10j) can be neglected, and we obtain n+(R) « 
C/R 2 , where C ~ 0.567 is the solution of the transcendental equation C = cxp(— C). 
Accordingly, the short-distance asymptote of e+(i?) reads 

(2-13) 

IraR^ mR z 

The 1 / i? 2 -potential does not have a characteristic length-scale and is at the origin of 
the Efimov effect. Note that in the case a = oo Eq. (|2.13l) holds for all R. 

2.4 Efimov effect in heavy-heavy-light system 

Let us assume that the light atom is in the symmetric state and that it follows adi- 
abatically the displacement of the heavy atoms (M/m ^> 1). Then, we can integrate 
out its motion and write down the Schrodingcr equation for the heavy atoms in the 
form 

[-V 2 n /M + e + (R)-E]<i>(R)=0. (2.14) 

At small R we can neglect E compared to other terms in Eq. (|2.14p and use the short- 
distance asymptote of the potential (|2.13[) . The radial Schrodingcr equation then reads 

{-d 2 /dR 2 + P/R 2 ) X (R) = 0, (2.15) 

where (3 = 1(1+1)— Q.lQM/m and x(i?) is the radial part of the Z-th angular momentum 
component of R<j)(H) . 

Equation (|2.15p is satisfied by the function R" ' , where the power v can take two 
values 

1/1,2 = 1/2 ± yJP + 1/4. (2.16) 

The case /? < —1/4 corresponds to the well-known problem of the fall of a particle 
to the center in an attractive R~ 2 potential (Landau and Lifshitz 1987). In this case 
a real solution of Eq. (|2.15[) is a linear superposition of R Vl and R V2 with complex 
conjugate coefficients. It can be expressed in the form 



X (R) oc VR cos[s \og(R/R )}, (2.17) 

where sq = — 1/4 and Rq is the three-body parameter, which fixes the phase in 

between the oscillations of R V1 and i?" 2 . 

The three-body parameter is sensitive to the short-range (high-energy) details of 
the interatomic potentials, i.e. two systems with different interparticle potentials are, 
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generally, characterized by different three-body parameters even if the corresponding 
scattering lengths are the same. The task of calculating Rq for a realistic potential is 
very challenging, but, most of the time, we need to know only a few general properties 
of this parameter. For example, for one and the same system it is reasonable to assume 
that the three-body phase is a continuous function of the potential depth, and in a 
sufficiently narrow vicinity of a two-body resonance it can be approximated by a 
constant. It then enters the zero-range theory as a fixed external parameter, whereas 
the scattering length changes in a very wide range. 

In the case of finite R e and a the wavefunction f|2 . 1 T[) remains a good approximation 
for distances R e < R < \a\. We see that it has approximately Nb ~ (sq/tt) log \a\/R e 
nodes, which is an estimate of the number of bound states (trimer states). This number 
becomes infinite in the limit a — > oo exactly as predicted by Efimov. 

Let us now look at the spectrum of trimers in the case a = oo. Equation (|2.14[) in 
this case reads 

(-d 2 /dR 2 + 13/ R 2 - ME)x(R) = 0. (2.18) 

Let us assume that there is a bound state with energy E = —Eq < 0. The wavefunction 
Xo(R) corresponding to this state is a linear superposition of ^fRJi So (iy/ M EqR) and its 
complex conjugate, where J is the Bessel function. The dependence of Eq on the three- 
body parameter is established by requiring that xo (R) be finite for large R ^> 1/ \/MEq 
and take the form (|2.17[) in the opposite limit. In principle, Eq itself can be used as a 
three-body parameter. 

Let us now consider a state given by the wavefunction xo(^R)- It obviously satisfies 
Eq. (|2.18|) with energy E = — X 2 Eq and decays at large R. The parameter A is not 
arbitrary - Eq. (|2.17p gives the condition sq log A = 7TO, where n is an integer. This 
is one of the consequences of the discrete scaling invariance. We understand now that 
all the Efimov trimers are self-similar and their spectrum is geometric: 

E n = E exp(27rn/s ). (2.19) 

Another manifestation of the discrete scaling invariance is the log-periodic depen- 
dence of three-body observables on the scattering length. This can be shown by rescal- 
ing the coordinates R = aR. In the new coordinates we have to solve the three-body 
problem for the unit scattering length and use the small- R three-body condition 

X (R) cx V^cob[« \og(R/R )}, (2.20) 

where the rescaled three-body parameter depends on the scattering length, Rq = Rq/ci. 
We then see that all the three-body observables should be log-periodic functions of 
a (besides the "trivial" dimensional scaling with a) with the discrete scaling factor 
exp(-7r/s )- 

As an example let us consider the s-wave atom-molecule scattering. In order to 
calculate the corresponding scattering length we take the / = component of Eq. (12.141) 
with E = — | eo | = —1/ma 2 , i.e. with the atom-molecule collision energy equal to zero. 
In the rescaled coordinates, R = R/a, this equation reads 



{-d 2 /dR 2 + Me+, a=1 (R) + M/m) X (R) = 0. 



(2.21) 
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Let us take the solution of Eq. (f2~2Tj) that at R < 1 behaves as R}l 2 + ia ° , it s long- 
range asymptote is a certain linear superposition of free solutions R and R 1 . We write 
it as A + BR, where A and B are complex coefficients, which could be in principle 
found numerically. Since Eq. (|2.21[) is real, the second linearly independent solution 
can be constructed as the complex conjugate of the first one. Then we should take their 
linear combination satisfying Eq. (|2.20[) and compare its long-range asymptote with 
the expression cx (R — a a d), where a a d is the rescaled atom-dimer scattering length. 
Performing this operation we find that 

ART is " + A* R is ° 

5ad(i?o) = ~ r I R* fr»o = * + C M*o log i?o + J) , (2.22) 

where £, Cj an d (5 are real numbers. We see that the nonrescaled atom-molecule scat- 
tering length has the form a a d = aaad(-Ro/ a ), i-e. it is indeed a log-periodic function 
of a (up to the dimensional prefactor). The resonances in the atom-dimer scattering 
length denote the emergence of a new trimer state each time the scattering length is 
multiplied by exp(7r/so). 

We emphasize once again that the discrete scaling invariance is a consequence of the 
fact that we assume Rq to be constant while we modify a. This is a good approximation 
close to a Feshbach resonance, however, theoretically we could rescale both the depth 
and the range of the potential so that Rq is proportional to a. Then, we would observe 
the usual continuous scaling, i.e. the self-similarity for any value of the scaling factor. 

2.5 Role of quantum statistics and masses of atoms 

In the Born-Oppenheimer approach we can also understand the role played by the 
quantum statistics of the heavy atoms. Indeed, the total three-body wavefunction is 
proportional to the product 0(R) i 0r ±(r). It should be symmetric (antisymmetric) 
with respect to the permutation R f± — R if the heavy particles are identical bosons 
(fermions) . Then the symmetry of </> depends on the choice of the light atom wavefunc- 
tion. Remember that ipn^(r) is symmetric and ipR,-( r ) is antisymmetric with respect 
to the permutation of the heavy particles. Since only the symmetric state leads to the 
induced attraction, the Efimov effect is expected for even angular momentum channels 
in the case of bosonic heavy atoms and for odd ones in the case of fermions. Of course, 
if the atoms are distinguishable, even if they are different hyperfine components of one 
and the same isotope, all angular momenta are allowed. 

Although the Born-Oppcnhcimcr approximation is supposed to work only for large 
mass ratios, exact calculations show that it captures quite well the essentials of the 
Efimov physics even for moderate M/m. In particular, the formula (see Section [2.41) 

s = sJO.lQM/m- (1 + 1/2) 2 (2.23) 

correctly predicts that the scaling factor A = cxp(7r/so) decreases with the mass ratio 
and increases with the angular momentum. The explanation of this phenomenon in the 
Born-Oppcnhcimcr picture is straightforward: an increase in the mass ratio strengthens 
the exchange l/i? 2 -attraction compared to the kinetic energy of the heavy atoms. In 
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M/m 

Fig. 2.2 Scaling factor exp(7r/so) (upper panel) and So (lower panel) versus M/m. The 
dashed lines are calculated by using the Born-Oppenheimer approximation and the solid 
lines are results of the exact approach of Sec. 13.21 Symbols show the scaling factors for 
several realistic three-atom systems. 

turn, the centrifugal barrier competes with the exchange attraction and prevents the 
heavy particles from falling to the center. 

In Fig. 12.21 we plot the scaling factor exp(7r/so) and sq versus M/m for I = and 
1 = 1. The result of Eq. (|2.23|) (dashed lines) is compared to the exact one (solid lines), 
the derivation of which will be discussed in Sec. 13.21 In the case of heavy bosons (I = 0) 
the Efimov effect always takes place. However, large mass ratios are more favorable for 
its observation as the scaling factor A rapidly decreases with M/m, which then allows 
for a larger number of log-periods within a fixed window of a. 

In the case of heavy fermions (I = 1) we see that the Efimov effect does not 
happen for M/m < 13.6 because the centrifugal barrier for the heavy fermions is 
stronger than the attraction induced by the light atom. The absence of the Efimov 
effect has advantages: First of all, the three-body problem in non-Efimovian cases 
is quite interesting on its own. For example, there can be trimer states which are 
qualitatively different from the Efimov trimers (Kartavtsev and Malykh 2007). Another 
striking feature of non-Efimovian systems is their collisional stability - the reason why 
a two-component strongly interacting fermionic mixture is long-lived even at unitarity, 
a = oo, whereas bosonic systems decay very rapidly near a Feshbach resonance. Let us 
now discuss the three-body recombination to deeply bound states - one of the most 
important loss mechanisms in a strongly interacting gas. 

2.6 Three-body recombination into deeply bound states: 
Efimovian case 

In the absence or in the case of a strong suppression of two-body inelastic channels 
a third body is needed for a pair of atoms to form a bound state. The binding en- 
ergy is then transfered into the kinetic energy of recombination products (an atom 
and a molecule). If the final bound state is deeply bound (typically of size R e ), the 
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products arc not ultracold and they escape from the trapped sample. The frequency 
of three-body recombination events in a unit volume is proportional to the number 
of all possible triples in this volume. For example, for a one-component Bose gas it 
is proportional to n 3 , where n is the density. The three-body recombination rate v 
is the frequency per atom, it is thus proportional to n 2 . The coefficient a roc in the 
equation v = a Icc n 2 is called the three-body recombination rate constant and has the 
dimension cm 6 s _1 . The loss rate equals 3^ as three atoms are lost at each event, and 
the density evolves according to the equation h = — 3im 3 . 

The recombination to deeply bound states can also happen in collisions of atoms 
and weakly bound molecules. In this case it is called collisional relaxation. Its frequency 
per unit volume is proportional to the product of the atomic and molecular densities, 
n a n m . Thus, the relaxation rate per atom equals a le in m , and per molecule, a xc \n a . 
The coefficient a re \ has dimensions cm 3 s _1 . 

Clearly, the task of calculating or, at least, estimating a rec and a ro i is very im- 
portant from the practical viewpoint. The knowledge of these constants allows one to 
tune experimental parameters and maximize the gas longevity. On the other hand, 
these parameters contain a lot of information on few-body correlations in the gas. In 
particular, the Efimov effect has so far been observed exclusively by measuring losses 
as a function of the scattering length. 

We now discuss the zero-range approach to this type of inelastic processes that 
happen at very short distances ~ R e . As an illustrative example let us consider the 
relaxation in atom-dimer collisions and calculate the corresponding rate constant a re i- 
It turns out that in the Efimovian case this problem can be solved rather elegantly by 
allowing the three-body parameter to be complex (Braatcn and Hammer 2007). Indeed, 
the solution of the heavy-atom Schrodingcr equation at distances R e <C R <C a is the 
sum of the incoming wave, y/R{R/R )-' s °, and the out going one, y/R~{R/ Ro) ls ° . In 
the presence of the relaxation to deeply bound states the outgoing flux should be a 
certain fraction of the incoming one. Theoretically this is modelled by introducing the 
so-called elasticity parameter 77* and writing Rq in the form Rq = \Rq\ exp(— t^/so). 
The modulus of Rq then fixes the relative phase of the incoming and outgoing waves, 
and rji, > ensures that the incoming flux is by the factor exp477* larger than the 
outgoing one. Everything that we have said about the three-body parameter holds 
also for the elasticity parameter. Namely, it is very difficult to calculate 77* ab initio, 
but in a sufficiently narrow vicinity of the resonance we can approximate it by a 
constant, which is inserted into the zero-range theory as a parameter. 

Now, let us consider a weakly bound molecule and an atom in a unit volume. The 
interaction energy shift for this system equals AE = 27ra a d//i, where jl w M/2 is the 
atom-molecule reduced mass and a a( j is the atom-dimer scattering length given by 
Eq. (j2.22j) . Once we introduce a finite elasticity parameter [Rq = \Rq\ exp(— ir/^/so) 
in Eq. (|2.22[) ] the atom-dimer scattering length acquires an imaginary part, which in 
turn leads to the shift of AE into the lower complex half-plane. We thus see that a 
single atom-dimer pair decays with the rate —2Im(AE), which in this case equals the 
relaxation rate constant 



47rlma ad 



47ra 



sinh 77* cosh 77* 



(2.24) 



Q^rel, bosons 



i 2 (s log \Ro/a\ + S) + sinh 2 77, 



cos 
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where £ and S are the same as in Eq. ()2.22j) . We see that a ro i, bosons/a is a log-periodic 
function of a. Quite naturally it reaches its maximum at the atom-dimer scattering 
resonances when the cosine in the denominator in the right hand side of Eq. (|2.24l) 
vanishes. 

In general, whenever we know the dependence of an arbitrary three-body observable 
on the real three-body parameter, we can analytically continue this function in the 
complex i?o-planc thus obtaining the result for the situation when the relaxation is 
present. As an exercise we propose to calculate the energies and lifetimes of the Efimov 
states in the case a = oo. 

The case 77* = 00 is quite peculiar. It corresponds to the complete absorption - 
only the incoming wave is present at short distances in the three-body wavefunction. 
In this case we recover the continuous scaling invariance, but three-body observables 
acquire a significant imaginary part. For example, the atom-dimer scattering length is 
complex and is simply proportional to a without log-periodicity. 

The real and imaginary parts of Rq can be found by fitting experimental loss data. 
Existing experimental measurements for various isotopes and Fcshbach resonances are 
consistent with values of rf in between 0.1 and 0.4. 

2.7 Three-body recombination into deeply bound states: 
non-Efimovian case 

The non-Efimovian case corresponds to (3 > —1/4 in Eq. (|2.16j) . Then the three-body 
wave function at distances R e <C R <C a is the sum 

X (R) oc A 1 (R/R e ) Vl + A 2 {R/R e ) v \ (2.25) 

where V\ > v-x- Expression (|2.25|) should be matched with the solution of the three- 
body problem at R < R e . To find the latter is challenging. However, in the absence of 
a three-body resonance the matching procedure implies that both terms on the right 
hand side of Eq. ()2.25|) are of the same order of magnitude, i.e. A\ ~ Ai- Then, at 
distances R ^> R e the second term can be neglected compared to the first one. Note 
that we use exactly the same argument when we neglect the interaction in between 
identical fermions in the ultracold limit: out of the two possible free solutions of the 
p-wave radial Schrodinger equation, R 1 and R~ 2 , we choose the former one (the one 
which increases with R). 

The zero-range approach to calculating the relaxation into deeply bound states in 
non-Efimovian cases is perturbative. It uses the unperturbed three-body wavefunction 
to predict the probability of finding three atoms at distances ~ R e and gives the 
functional dependence of a re i on the scattering length. This means that if the relaxation 
rate constant is known (measured) for a certain a, one can predict its value for any 
other a^> R e . 

As a specific example let us estimate the atom-dimer relaxation rate in p-wave colli- 
sions of heavy fermions with heavy-light molecules in the Born-Oppenheimcr approxi- 
mation. In this case the light atom is in the state i/sR. t +(r), which means that at R <C a it 
is always close to the heavy fermions. Then we just need to calculate the probability to 
find the heavy atoms at distances R ~ R e . In order to do this we take the short-distance 
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(R < a) asymptote of the heavy-atom wavefunction, ^(R) ps CR Ul ~ l cos 6>k, r, and 
match it with the long-distance one, $(R) ps V2 sin(kR) . The latter is properly nor- 
malized to a single atom-molecule pair in a unit volume (we assume that we are 
far from any atom- molecule p-wave resonance), k is the relative atom-dimer mo- 
mentum, and #k,R is the angle between k and R. Matching the short- and long- 
distance asymptotes at R ~ a we obtain C ~ ka 2 ~" 1 , which gives the probability 
Pr<r ~ (R e / 'a) 2vi ~ 2 '(ka) 2 In order to get a dimensional estimate of a re i,fermions 
we can multiply Pr<r c by the frequency of recombination processes that would take 
place if all three atoms are confined to distances ~ R e . This frequency is of order 
1/mR 2 . We thus obtain 

arci,fcrmio„s = (R e / m)(R e / a) 2 ^ 2 (ka) 2 . (2.26) 

The three factors in Eq. (|2.26|) are interpreted as follows: The factor R e /m is of 
the order of the relaxation rate constant for s-wave collisions of atoms and deeply 
bound molecules (with the size ~ R e ). The factor (R e /a) 2l/1 ~ 2 is the suppression factor 
that comes from the fact that the atoms have to tunnel under the effective repulsive 
potential, which is the sum of the centrifugal barrier and the exchange attraction. The 
power v\ depends on the mass ratio. The factor (pa) 2 corresponds to the low-energy 
Wigner law for reactions with unit angular momentum (I — 1 in this case). 

Accordingly, the relaxation rate constant for s-wave atom-molecule collisions has 
the form (|2.26[) without the last factor. Unfortunately, in this case the Born-Oppenheimer 
approximation does not work because the state ipii,-( r ) is unbound at R < a. The 
problem should then be solved without relying on the adiabatic approximation (see 
Sec. 133)1 . 

2.8 Role of resonance width 

As we have explained in Sec. 12.41 the Efimov effect is related to the effective 1/R 2 - 
attraction emerging in a three-body system at distances R < \a\. Let us now discuss 
how this effective potential is modified in the case of a narrow resonance. First, consider 
the regime of small detuning, R e -c R* <C \a\. In this regime the right hand side of 
Eq. (|2.10[) can be neglected at distances R* -C R -C |a|, and we recover the 1 /Mis- 
behavior of the effective potential (|2.f 3[) . However, at distances R e -C R <C R* the 
parameter k + (R) is determined mostly by the last terms on each side of Eq. (|2.10|) . 
We thus find that k 2 + (R) ps 1/R*R and 

e + (R) ps -l/2mR*R, R e < R < R* < |o|. (2.27) 

The Coulomb potential (|2.27|) is qualitatively different from the inverse-square one 
(|2.13[) . In particular, in this potential heavy atoms do not fall to the center. Moreover, 
close to the origin, R -c (m/M)R* , the heavy-atom wavefunction behaves similarly to 
the noninteracting case, i.e. at these distances the kinetic energy operator in the heavy- 
atom Schrddinger equation is dominant. Therefore, we do not have to worry about the 
boundary condition at the origin (as in Sec. 12. 71 we just choose the solution that grows 
faster). The wavefunction is thus uniquely defined. Then at distances R ~ R* it can 
be matched with the Efimov-likc wavefunction (|2.17j) , the three-body parameter being 
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determined by R* . Remarkably, in this case three-body observables depend only on 
the two-body parameters a and R* . 

It is instructive to discuss more qualitatively what happens with the system as R 
decreases from R^> R* to R <C R* . One can see that in the former case the effective 
energy-dependent scattering length a(E) introduced in Eq. (|1.18|) is larger than R 
and, therefore, drops out of the problem as in the wide-resonance case with a = oo. 
However, for R ~ R* the energy of the light atom becomes sufficiently detuned from 
the resonance, so that a(E) is comparable to i?, and the effective potential acquires 
a characteristic lengthscale and is no longer that steep. Another important point to 
mention is a qualitative change of the wavefunction at R ~ R*. By using Eqs. (|1.27|) 
and (|2.6[) it is straightforward to show that for R 3> R* the light atom occupies 
predominantly the open channel, whereas at distances R <C R* the open channel 
occupation is negligible. Therefore, at these distances the wavefunction of the system 
to the leading order describes the free motion of a heavy atom and a closed channel 
molecule. A small open channel occupation can be treated perturbativcly. 

In the regime of intermediate detuning, R e <C \a\ <§C R* , the effective potential is 
nowhere proportional to 1/R 2 . At distances R < \a\ it is approximated by the Coulomb 
potential (|2.27[) . Then in the case a < there are no bound states for the light atom if 
R > \a\ and the potential e+(R) terminates at the point R = \a\ where it reaches the 
three-atom continuum. For positive a, the effective potential changes its shape from 
the Coulomb one for R < a to the constant e+(R) ~ — |eo| for R > a. Accordingly, in 
both cases trimer states appear only when the characteristic Bohr radius, oc (m/M)R* , 
corresponding to the potential (|2.27[) is smaller than \a\. 

The last point that we would like to touch upon in this section is inelastic losses 
in few-body systems near a narrow resonance R* . The approach to this problem is 
similar to the one discussed in Sec. 12.71 for non-Efimovian systems. In fact, in the 
narrow-resonance case the probability of relaxation or recombination to deeply bound 
states can always be related to the relaxation rate constant for collisions of atoms and 
closed channel molecules. Indeed, the relaxation is a local process, it requires three 
atoms to approach each other to distances ~ R e -C R* and, as we have argued above, 
at these distances the wavefunction describes the relative motion of an atom and a 
closed channel molecule. We will return to this question in Sec. 13.51 More details on 
the few-body problem near a narrow resonance can be found in (Petrov 2004; Wang 
et ai. 2011a; Levinsen and Petrov 2011). 
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The method of Skorniakov and 
Ter-Martirosian (STM) for few-body 
problems with resonant short-range 
interactions 



In spite of numerous advantages and simplifications provided by the Born-Oppcnheimer 
approximation we still have to discuss methods suitable to quantitatively describe the 
case of comparable masses and even some highly mass-imbalanced ones (for example, 
the s-wave scattering of a heavy fermion and a weakly bound heavy-light molecule 
when the symmetry forces the light atom to be in the state "0R,- ( r ) which does not 
exist at small R). 

A natural generalization of the Born-Oppenheimer approximation is the adiabatic 
hyperspherical method, in which, after separating out the center-of-mass motion, one 
introduces the coordinate system of hyperradius and hyperangles, {p, Cl}. The hyper- 
radius p is the square root of the sum of squares of all interparticle distances with 
mass-dependent weights, and all other coordinates of the system can be written as a 
set of dimensionless hyperangles f). The adiabatic idea in this case is to consider p 
formally as a slow coordinate and Cl as fast ones. Accordingly, one fixes p and diago- 
nalizes the hyperangular part of the Schrodinger equation thus obtaining the channel 
potentials. The second step is then to solve the hyperradial part, which is a set of 
second-order ordinary differential equations for each channel. The difference from the 
simple Born-Oppenheimer approach is that all these differential equations are now 
coupled by nonadiabatic matrix elements. In practice one truncates the set of these 
coupled equations according to a given accuracy goal. Usually only a few channels 
suffice to obtain convergent results. 

An advantage of the adiabatic hyperspherical approach is that it is quite general. 
It can be used for finite- and long-range interaction potentials, one can easily include 
three-body forces, etc. In fact, the complexity of the method is equivalent to the 
complexity of the initial few-body Schrodinger equation and is practically independent 
of the form of the interatomic potentials. Let us consider an .ZV-body problem in three 
dimensions (TV > 3). Then we can go to the center-of-mass reference frame and separate 
out the three global Euler angles by using the rotational invariance. We end up with 
3iV — 6 degrees of freedom, 3N — 7 of which are hyperangular. Thus, for N = 3, the 
hyperangular space is two-dimensional, and each additional particle adds 3 degrees 
of freedom, which makes the hyperangular part of the calculation quite challenging. 
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For the details on the adiabatic hyperspherical approach see, for example, (Lin 1995; 
Nielsen et al. 2001) 

Another approach, appropriate for systems with short-range interactions, was first 
introduced by Skorniakov and Ter-Martirosian (Skorniakov and Ter-Martirosian 1957). 
For N = 3 it leads to a one-dimensional integral equation in momentum space (the 
STM equation), and its generalization to TV > 3 gives an integral equation for a 
function of 3iV — 9 coordinates. The STM equation can be obtained by using the 
Effective Field Theory (EFT) (Bcdaque et al. 1999) and other diagrammatic techniques 
(Brodsky et al. 2005; Levinsen and Gurarie 2006). We should mention that the famous 
Faddeev-Yakubovsky few-body method can be understood as a generalization of the 
STM approach to the case of finite-range (and long-range) potentials (Faddeev 1960; 
Yakubovskii 1967; Faddeev and Merkuriev 1993). 

We will now derive the STM equations in coordinate space directly from the 
Schrodinger equation. This derivation clearly shows how one can integrate out all 
short-range physics and use the fact that mostly everywhere in space the motion of 
the atoms is free. 

3.1 STM equation 

Let us first demonstrate the general idea behind this method. Assume that we are 
in two-dimensional space and we solve the problem of scattering by a curve S (see 
Fig. 13. ip . Namely, we solve the free Schrodinger (or Hclmholtz) equation 

(-V^-£)V = (3.1) 

with the boundary condition 

[di/>/dn]/ip = F(x), (3.2) 

where x is the natural (by the arc length) parametrization of S, F is a known function, 
and [dip /dn] is the sum of the normal outward derivatives at the point x of the curve. 

An efficient way of solving this scattering problem is to introduce an auxiliary 
function f(x) defined on the boundary and look for the solution in the form 

Hp) = MP) + [ G E [P - Ps (x)]f(x)dx, (3.3) 
Js 

where Ge(p) is the Green's function of Eq. (|3.1|) . Ps{%) is the coordinate on the curve, 
and ipo(p) is the incoming wave, which is a free solution of Eq. (13. ip . 

Clearly, the function (|3.3p satisfies Eq. (|3.ip . We just have to ensure the boundary 
condition (|3.2p . It can be shown (by choosing a proper contour around the curve 
clement dx and applying the Gauss-Ostrogradsky theorem) that [dip/dn] = f(x). 
Then Eq. (|3.2p gives the equation 

MPs(x)} + f G E [ Ps (x) - p s {x')]f{x')dx' = 1/F(x). (3.4) 
Js 

Note that Eq. (|3 .4[) is one-dimensional in contrast to the original two-dimensional 
Schrodinger equation ()3. 1|) . In fact, it is much more suitable for analytical and espe- 
cially numerical calculations as this reduction of the configurational space allows for 
rapid computational schemes. 



STM equation 23 





Fig. 3.1 Two-dimensional scattering of a particle by a curve. 

The idea of using free Green functions is quite natural when one deals with bound- 
ary value problems. It can be encountered in electrostatics, hydrodynamics, problems 
of heat diffusion, etc. We can mention the so-called method of boundary elements used 
for calculating spectra of quantum billiards (Berry 1981; Bohigas et el. 1984) and based 
on the Korringa-Kohn-Rostoker method (Korringa 1947; Khon and Rostokcr 1954) in 
solid state theory. Let us now discuss how one can do similar things for the few-atom 
problem. 

Consider three atoms with coordinates r* and masses rtii, where i = 1,2,3. We 
separate out the center of mass motion and introduce the three sets of rescaled Jacobi 
coordinates (see Fig. 13. 2[) : 

(mjTj + m k Yk)/{m.j + m k )], 

-r 3 -), (3.5) 

where k} are cyclic permutations of {1, 2, 3}, i.e. {1, 2, 3}, {2, 3, 1}, and {3, 1, 2}, 
\H = mjmk/(m,j + m^) is the reduced mass for atoms j and k, and p,i = mi(mj + 
nik) / (rrii + nij + rrik) is the reduced mass for the relative motion of atom i with respect 
to pair jk. The purpose of the mass-dependent rescaling is to express the kinetic energy 
operator in the symmetric form — V£ 4 — V yi . 

Let us choose to work in the first coordinate system (i = 1 in Eqs. p.5[l ). It is 
related to the other two by the equations: 




Xl = Xi^ 2 (x 2 ,y2) = -(V^lA i 2/TO3)x 2 + VM2/Aiy2, 

yi = Yi^ 2 (x 2 ,y2) = - V M1/A2X2 - (VMiM2/»i3)y2, (3.6) 

and 



xi = Xi^ 3 (x 3 ,y 3 ) = -( v / /"iA i 3/?«2)x3 - VM3/Aiy3, 

yi = Yi < _ 3 (x3,y 3 ) = v//ii/// 3 x 3 - ( v / /"iA i 3/m 2 )y3. (3.7) 
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Fig. 3.2 Three sets of coordinates for the three-body problem. 



Now the original three-body problem can be represented as a single particle scat- 
tering in the six-dimensional space by a potential which is localized on the three- 
dimensional hyperplanes defined by the equations y$ = 0, i = 1,2,3. Everywhere 
outside these hyperplanes the three-body wavefunction satisfies the free-motion six- 
dimensional Schrodinger equation 

(-V» 1 -Vj 1 -£?)*(xi I y 1 )=0. (3.8) 

As in the usual scattering theory we write ^(xi, yi) = vE^xi, yi)+vE , S cat(xi, yi), where 
( x i , y l ) is the incoming wave - a free solution of Eq. (|3.8j) without singularities (as if 
there were no interactions), and 4 , sca t(xi, yi) is the scattered wave which contains 1/y- 
singularities at each of the boundaries. We now introduce auxiliary functions defined 
at the boundaries, /i(xj), and write ^(x^yi) in the form 

*(xi,yi)=* (xi,y 1 )+ / G B ( v /(x 1 -x') 2 +y?)/i(x')d 3 a ;' 

+ E / °e(VWi -Xi^(x',0)] 2 + [yi - Yi^ i ( x ',0)] 2 )]/ i (x')dV, (3.9) 

i=2,3 J 

where Ge is the Green's function of Eq. (|3. 81) : 

i 8tt 3 X 2 ' " < U 1 

jgW^) p > n E ^ Q 4n x 
— le^x* — > ^ > u 



By construction, the wavefunction (|3.9p satisfies the Schrodinger equation (|3.8[) 
for arbitrary boundary functions /j(xj). As in the previous example this freedom is 
removed by boundary conditions. We first discuss the wide resonance case. Having in 
mind the Bethe-Peierls boundary condition (|1.9[) let us study the singular and regular 
terms of ^(xi,yi) close to the boundary y\ —¥ 0. In this limit ^o( x i,yi) and the two 
terms in the sum in the second line of Eq. (|3.9[) are not singular. We can thus set y± = 
there. In the remaining term we make the 1 / yi-singularity explicit by subtracting and 
adding the quantity 



/i(xi) f G B ( v /(x 1 -x') 2 + y 2 )dV = /ifxOe-^^^i, (3.11) 
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where, for E > we use the convention y/—E = —i\f~E. We thus get the following 
asymptotic expression for ^(x^yi): 

*(x ljyi ) -> (l/4 7 r)[(l/y 1 -y=^)/ 1 (x 1 ) + 47rf (xi,0)- J L 1 , E {/}(x 1 )] J (3.12) 

Hi->0 

where we have expanded the right hand side of Eq. (|3.11j) to the next-to-leading order 
at small y\, {/} = {fi,h, /a}, and Li t E is the integral operator defined by 

£i,b{/}(xi) = 4tt /" G B (|xi - x'DI/^xx) - /i(x')]dV 

-4tt ^ /*G £ [^/a;5+a;' 2 +2^iM(l-w/mi)xix']/ i (x')dV. (3.13) 

1=2,3"' 

Comparing Eqs. (|3.12[) and (|1.9|) we arrive at the integral equation 

Li.B{/}(x 1 ) + (V^B-l/v / 2MTai)/i(xi) = 47r* (x 1 ,0), (3.14) 

where ai is the scattering length corresponding to the interaction of atoms 2 and 3. 

Applying the same procedure to the boundaries y 2 = and y 3 = we can de- 
rive two additional equations, the left hand sides of which can be obtained from 
that of Eq. (|3.14p by cyclically permuting subscripts 1,2, and 3. The right hand 
sides of the new equations read, respectively, 47r\f , o[Xn_2(x2, 0), Yi < _ 2 (x2, 0)] and 
47T<I'o[Xi^3(x3, 0), Yi < _3(x3, 0)]. We thus get three coupled integral equations for de- 
termining the boundary functions fi and we have substantially reduced the configura- 
tional space of the problem. In fact, the operators Li t E conserve angular momentum, 
and we can expand fi in spherical harmonics to work only with a set of one-dimensional 
integral equations. Solution of these equations then gives us the wavefunction ^ by 
virtue of Eq. (|3.9[) . Vice versa, if we know VP, we can always obtain fa by looking at 
the coefficients in front of the 1 /^-singularities of \J> at small j/j. Namely, 

fifa) = 4tt lim w *. (3.15) 

2/i->0 

Thus, {/} contains the same information about the system as the wavefunction ^. 
This is one of the key advantages of the zero-range approximation. 

It is sometimes useful to go to momentum representation introducing the Fourier 
transform /fc(p) = / /fe(x) exp(— ipx)d 3 x. Then the operator Li j _e{/} reads 




LiMf}(pi) = W-e+pI - v^)/i( P i) 

/.(pW 



i=2,3 



p' 2 +Pi + 2y/n 1 /m(l - fMi/mi)pip' - (jJLi/fll)E 



(3.16) 



Skorniakov and Ter-Martirosian derived Eq. (|3 . 14[) in momentum representation in the 
case of equal masses in order to solve the problem of neutron-deutron scattering (Sko- 
rniakov and Ter-Martirosian 1957). We will refer to Eq. (|3.14j) as the STM equation. 
Let us now illustrate the power of the method by considering some concrete examples. 
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3.2 The Efimov effect and determination of s Q 

We return to the problem of Sec. 12.51 and show how to determine the parameter so 
without relying on the Born-Oppenheimer approximation. The system we are inter- 
ested in consists of two identical atoms (bosons or fermions) of mass mi = m 2 = M 
interacting resonantly (scattering length a) with another atom of mass m 3 = m. Let 
us neglect the (nonresonant) interaction between the identical atoms. Then \£ is not 
singular for y 3 — > and therefore / 3 = 0. By symmetry the boundary functions /i 
and / 2 are equal to each other (bosons) or have different signs (fermions). Indeed, 
directly from Eq. Q3.5P we see that the permutation ri <=^ r 2 is equivalent to Xi T± x 2 
and yi T± — y 2 . Thus *(xi,yi) = ±\&(x2,— y 2 ) and by virtue of Eq. (|3.15|) we have 
/(x) := /i(x) = ±/ 2 (x), where the upper sign stands for bosons and lower - for 
fermions. Clearly, should have the same property. Finally, the three equations of 
type p.!4[) are reduced in this case to the single one, 

Wto + - l/v / 2Ma)/(x) = 47r*o(x, 0), (3.17) 

where 

L E /(x) = 4irJ {G B (|x - x'|)[/(x) - /(x')] T G E (^x 2 + x> 2 + 2xx' sin0)/(x')}dV, 

(3.18) 

/i = hi = /i 2 = mM/(m + M), and sin</> = M/(m + M). 

At small x <§; {a, we can neglect all terms in Eq. ()3.17|) except the integral 

one and in the latter we can use the zero-energy asymptote of the Green's function 
(|3.10|) . This means that the small- a; asymptote of /(x) should satisfy 

L s =o/(x) = 0. (3.19) 

We first consider spherically symmetric /(x) = f(x). In this case averaging over the 
angles of x' (we integrate over the solid angle and divide by Air) we obtain the zero 
momentum component of the operator Le=o : 



4 
ir 



t' 2 Ht' f3 20) 

(x 2 - x' 2 ) 2 ^ (x 2 + x' 2 ) 2 - 4x 2 x' 2 sin 2 ■ ' 

The first integral in Eq. (j3~20|) [and also in Eqs. (j3~T3| and ([3"T8]) ] is taken in the 
principal value sense (Petrov et al. 2005a). We see that the operator (|3.20[) is scaleless 
and has the property 

L B=0 j= x v = AfcoM^- 1 . (3.21) 

We thus look for the solution of Eq. (|3.19|) in the form f(x) = x v , where v is a root of 
A/ = o( i/ )- The region of convergence of the integral on the left hand side of Eq. (|3.21l) is 
—3 < Re(^) < 1 in the case of bosons and —5 < Rc(is) < 3 for fermions. The function 
\i=o(v) is given by 

iris 2 sin[fl!)(^ + 1)1 
\l=o(v) = —(v + 1) tan — =p , ^ ' 3.22 
2 sm(20) cos(7r^/2) 

The integrals of type Le=o.i=qx u can be taken by using the complex analysis. We cut 
the complex x-plane along the positive real axis and choose the integration contour 
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embracing this cut. Then the contour can be blown to infinity, and the integral is 
determined by the pole residues of the integrand. 




v 



Fig. 3.3 \i =0 (v) versus v in the case of two identical bosons (solid lines) or two identical 
fermions (dashed lines) interacting resonantly with a third atom of mass m. The Efimov 
effect takes place only for bosons. The star and cross denote the value v\ ~ f .f6622 for equal 
masses and v\ ~ 2.0f93 for the case of two 40 K and one 6 Li atoms. 

In Fig. 13. 31 we plot \i = q(v) in the bosonic (solid lines) and fermionic (dashed lines) 
cases. In the fermionic case for any mass ratio the equation X(y) — has two solutions 
v\ > vi in the region of convergence mentioned above. These roots are real which 
means that we are in the non-Efimovian case discussed in Sec. 12.71 Accordingly, in 
the absence of a three-body resonance f(x) oc x" 1 . Note that we have just solved the 
problem that could not be solved in the Born-Oppenheimer approximation even at 
very large mass ratios. Indeed, in the case of heavy identical fermions the symmetry 
I = corresponds to the light atom in the state ip- , which is unbound at these small 
distances. 

The case of identical bosons in the Born-Oppenheimer language corresponds to the 
wavefunction ip + for the light atom and zero total angular momentum (see Scc. l2.5[) . In 
this case the Born-Oppenheimer approach predicts the Efimov effect. Now, studying 
A; = o (v) one can show that the effect takes place for any mass ratio as the roots of this 
function are complex conjugate, v\ 2 = — 1 ± is$, where sq is real. Accordingly, the 
three-body parameter is necessary in order to fix the ratio in between the coefficients 
in the linear superposition of x Vl and x" 2 . The parameter sq plotted in Fig. 12.21 (solid 
line) is found by solving the equation A; = o(— 1 + iso) = 0. The increase of so with the 
mass ratio is consistent with the Born-Oppenheimer picture. 

As a historical remark we note that the neutron-deutron s-wave scattering problem 
considered by Skorniakov and Ter-Martirosian for total spin 7 = 3/2 corresponds to the 
fermionic case discussed above and for total spin I = 1/2 - to the bosonic case. Indeed, 
in the latter case the neutron spins are antiparallel and their orbital wavefunction is 
symmetric as for bosons. It was later pointed out by Danilov (Danilov 1961) that for 
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I = 1/2 the solution of the STM equation p,17|) is not unique and an additional 
parameter is necessary. 

We should also note that the bosonic case with equal masses is characterized by 
relatively small sq « 0.414 leading to very large scaling factor exp(7r/so) ~ 1986.1. This 
fact, although not disturbing for theorists, causes apparent practical difficulties. The 
reason for this "weak" manifestation of the Efimov physics is that there are only two 
resonant interactions out of three. However, as we have already mentioned in Sec 12.51 
the scaling factor rapidly decreases with the mass ratio (see upper panel of Fig. 12.21) . 
which makes highly mass-imbalanced heteronuclear mixtures practically valuable for 
studies of the Efimov effect and discrete scaling invariance. 

The case I = 1 is treated in the same manner as I = 0. We obtain the operator 
Le=o,1=i by integrating out the angular dependence of /(x) = 3/(x)xn in £e=o/(x) 
(we multiply it by xn, integrate over the solid angle, and divide by 47r): 



L E =o,i=if{x) 
/(*') 



1 

7T 



4x' 2 [f(x)-f(x')} f{x') 



[x r 



2xx' 



T 



2x 2 sin 



4xx'(x + x ) sin0 
+ x' 4 + 2x 2 x' 2 cos(20) 



+ log ■ 



(x + x') 2 

2xx' sin 



log 



x 



X 



12 

X 



2xx' sin < 



x — x 
dx'. 



(3.23) 



The operator (|3.23|) has the same scaling property p. 211) with 



, v(v + 2) 7w usm((f)) cosUu + 1)0] — sin(i/0) 

A; — i (v) = cot — =F ^ ■ 

W v + 1 2 T (u + l)sin 2 (0)cos((/))sin(7r^/2) 



(3.24) 



The region of convergence of the integrals in Eq. (|3.23|) is —4 < Rc(^) < 2 for bosonic 
and for fcrmionic symmetries. 

In Fig. 13.41 we plot the function \i = \(v) for bosons and for fermions. In the Born- 
Oppenheimer language this symmetry (1 = 1) corresponds to the antisymmetric heavy- 
atom wavefunction. In the bosonic case, besides leading to the centrifugal barrier, it 
also forces the light atom to be in state tp_ (leading to an effective repulsion for the 
heavy atoms). This explains the absence of the Efimov effect for bosons in this case. 
For fermions we have the competition of the centrifugal barrier and the attraction due 
to the exchange of the light atom in state For mass ratios M/m > 13.6 the latter 
wins and the Efimov effect manifests itself in the complex roots of A/ = i (y ), which have 
the same structure as in the case of zero angular momentum, = —1 ± isq. The 
quantity sq for the I = 1 case is also plotted in Fig. 12.21 along with the result of the 
Born-Oppenheimer approximation. 



3.3 Analytical approach at zero total energy 

In this section we would like to present an approach which allows one to calculate all 
zero-energy three-body obscrvables analytically. The approach was introduced for the 
fermionic non-Efimovian case in (Petrov 2003) and has been recently generalized to 
Efimovian cases in heteronuclear mixtures (Helfrich et al. 2010). Here, we apply it to 
three identical bosons. Alternative analytical methods of treating this system can be 
found in (Macek et al. 2005; Macek and Ovchinnikov 2006; Gogolin ct al. 2008; Mora 
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Fig. 3.4 Xi-i(v) versus v in the case of two identical bosons (solid lines) or two identical 
fermions (dashed lines) of mass M interacting resonantly with a third atom of mass m. For 
this symmetry (I = 1) the Efimov effect is absent in the bosonic case. For fermions it takes 
place only for mass ratios larger than the critical one, (M/m) c « 13.6. 



et al. 2011a). We should mention that our current understanding of the three-boson 
problem is strongly influenced by the papers of Nielsen and Macek (1999), Esry et 
al. (1999), Bedaque et al. (2000), and Braaten and Hammer (2001), who analyzed 
the process of three-body recombination in this system numerically and arrived at a 
number of important conclusions. 

Close to a Fcshbach resonance a Bose gas suffers from three-body recombination - 
formation of weakly (for large positive a) or deeply bound molecules, the binding en- 
ergy being transfered to the kinetic energy of the products - a molecule and remaining 
atom. Let a s , ad{a > 0), and ctd{a < 0) stand for the rate constants for recombina- 
tion to a weakly bound state and deeply bound states on the positive and negative 
sides of the resonance, respectively. In the low temperature limit, as long as the de 
Broglie wavelengths are larger than |o|, these constants are practically temperature 
independent and can be written as products of ha /m and dimensionless log-periodic 
functions of |a|/i?o [different for a s , ad(a > 0), and ctd{a < 0)]. For fixing the phase 
of these log-periodic functions the following practical notation is used in literature: 
the value of a > where a s /a reaches its minimum is denoted by a*o [defined mod- 
ulo exp(7r/so)]. Another experimentally relevant reference point is the maximum of 
otd(a, < 0)/a 4 , which is by definition reached at a = a_ < 0. In fact, this is the 
point where a Efimov state crosses the three-atom threshold. We will show below that 
|a_|/o»o = exp(7r/2s ). 

Let us first discuss the case a > and temporarily adopt the units h = m = a = 1. 
For the problem of three-body recombination we now consider Eq. (|3.14j) and choose 
the incoming wave \l/o as a symmetrized superposition of plain waves normalized to a 
volume V. In the region of space relevant for recombination it can be approximated 

by 
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*o = V6V~ 3/2 , (3.25) 

where the factor a/6 is due to the fact that we consider cold bosons but not in the same 
quantum state (not condensed). Let us mention that Eq. Q3.14p with E = can also 
be used for the problem of atom-dimcr scattering just below the breakup threshold 
(E = — 0). In this case the atoms cannot move freely at large distances and 'J'o should 
be set to zero. 

Using the bosonic symmetry and Eq. (|3.15|) it is straightforward to show that 
/(x) := /i(x) = /2(x) = /s(x) (we work in the first coordinate system, {xi,yi}, 
and omit the subscript). Equation (|3.14[) is thus a single three-dimensional integral 
equation. Moreover, we can set /(x) = f(x) consistent with the angular independence 
of and with the fact that all processes with higher angular momenta are suppressed 
at low temperatures. The resulting STM equation reads: 

(£e=o,i=o - l)/(a:) = 4tt* , (3.26) 

where the integral operator Le=o,i=i is obtained in the same manner as (|3.20[) . Omit- 
ting the subscripts, we get 



4 f° 
Lf(x) = - \ 
* Jo 



f(x) - fix') 2f(x') 



(x 2 — x' 2 ) 2 (x 2 + x' 2 ) 2 — x 2 x' 2 ^ 

Let us now discuss the structure of possible solutions of Eq. f|3 . 26|) . Obviously, f{x) 
is a sum of a particular solution of the inhomogeneous equation p.26|) and a general 
solution of the homogeneous equation 

(L - l)x(x) = 0. (3.28) 

Physically, Eq. (|3.28[) describes the atom-dimer channel just below the dimer breakup 
threshold (Wo = 0). Therefore, at distances x ^> 1 the function xi x ) 1S a linear com- 
bination of exp(ix)/x and exp(— ix)/x. Indeed, consider the atom-dimer wavefunction 
ty(x,y) = 4>b{y) exp(±Ax)/x. Substituting it into Eq. (|3.15|) and using the correctly 
normalized dimer wavefunction 



My) = cxp(-y)/V27rj/ J (l/y - 1)/V2vr (3.29) 

we find that the corresponding contribution to x( x ) equals (87T) 1 / 2 exp(ix)/x. 

Our aim now is to solve Eq. (|3.26[) , separate the large- a; asymptote / cx e~xp(ix)/x, 
and relate the coefficient in front of it to the three-body recombination rate constant. 
This problem can be solved analytically by using the property ()3.21[) of the operator 
L. In the case of three identical bosons the integral converges in the region —3 < 
Rc(^) < 1 and the function A; = o(^) is given by (we again omit the subscript) 

, . . . , iris 8 sin[7r(f + l)/6l 

X(v) = — (i> + 1) tan — = — LA Jl± m 3.30 

w v ' 2 V3 cos(tt^/2) v ; 

The roots of (|3.30p are complex conjugate, v\ 2 = — 1 ± isq, where sq satisfies 

s cosh(7rs /2) - 8 sinh(7rso/6)/V3 = 0. (3.31) 

The solution is so ~ 1.00624 which leads to the famous scaling factor exp(7r/so) ~ 22.7. 
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At short distances the operator L in Eq. (|3.28[) dominates over 1, and any solution 
of this equation should be a linear superposition of \ oc x +iS ° and its complex 
conjugate. From now on we will use the notation \ for the solution of Eq. (|3.28[) with 
the following asymptotes: 



X-( x ) — { x -l e ix+i<r-h _|_ x -l e -ix-ia+h x ^ I 



where A is a complex number, and er and h are real numbers. The physical solution 
of Eq. (|3.28p . i.e. the one corresponding to a given three-body parameter, is expressed 
as 

Xe(x)=e ie x(x)+e- ie X *(x), (3.33) 

where we have introduced the three-body parameter 9 (a complex number with imag- 
inary part T]*). 

The normalization in Eq. p.32[) is chosen such that 

/•OO 

(px(p x )\p'x(p' x )) = / Px{px)p'x{p'x)x 2 dx = 2tt6(p - p'). (3.34) 



The first equality in Eq. p.34[) is our definition of the scalar product (note the absence 
of the complex conjugation), and the second equality follows from the fact that px(px) 
and p'x^P' x ) are eigenfunctions of the symmetric operator L corresponding to the 
eigenvalues p and p' . They are orthogonal for p ^= p' and their scalar product in the 
vicinity of p = p 1 can be worked out in the same way as in (Landau and Lifshitz 
1987, sec §21). A simple change of the integration variable in Eq. (|3.34[) leads to the 
completeness condition 

poo 

p 2 x{px)x{P x ')dp = 2tt8(x — x')/x 2 . (3.35) 

Equations (|3.34j) and (|3.35[) allow us to construct the integral operator (L — l) -1 
needed to solve Eq. (|3.26[) . In order to avoid problems with divergence of the corre- 
sponding integrals let us introduce an auxiliary function go{x) related to /(x) by 

/(*) = 4tt*o[-1 - A(0)/.t + A(0)A(-l)so(x)]. (3.36) 

Substituting this expression into Eq. (|3.26|) and using Eq. (|3.21[) we find that go(x) 
satisfies the equation (L — l)go(x) = x~ 2 . Applying the operator (L — 1) _1 to x~ 2 wc 
obtain the following particular solution 



9o(x) = ^ / 



where the first integral is defined as 



x(y)ydy 2mxx(x) 



y — x — iO 1 — exp(— 2ttsq) 



(3.37) 



x{z)dz = lim / x(z)z c rfz. (3.38) 

e— >+0 



The rule of going around the pole in the second integral and the numerical coefficient 
in front of the second term in the square brackets on the right hand side of Eq. (|3.37|) 



32 The method of Skorniakov and Ter-Martirosian (STM) for few-body problems with resonant short-range interactions 

regulate the entry of x(x), which can be arbitrary, into the particular solution go(x). 
Using this freedom we choose these parameters in such a way that go (x) does not con- 
tain oscillating terms proportional to x~ 1+ls ° at small x. Direct calculation shows that 
in the limit x — >• the right hand side of Eq. (|3.37[) equals go(x) ~ [J Q °° x(z)dz] 2 /2ttx 
to the leading order in x. On the other hand, according to Eq. (|3.21[) the same quantity 
in the same limit can be written as go(x) = (L — l) _1 a; _2 w L _1 x~ 2 = 1/[A(— l)x], 
which leads to the result 



Another consequence of our choice of the particular solution p.37[) is that removing 
the oscillating terms from go{x) makes it real, since any imaginary part of go would 
necessarily be a solution of the homogeneous Eq. (|3.28[) . Therefore, go would have 
oscillations at short x, the absence of which we have ensured. Clearly, the function / 
obtained by virtue of Eq. p.36[) is also real. Moreover, the property (|3.21[) ensures that 
/ = o(l) at small x, i.e. its Taylor expansion starts with x 1 , at least. Therefore, this 
solution of Eq. (|3.26l) is not sensitive to the short-range physics and does not depend 
on the three-body parameter. 

Integrating Eq. (|3.37[) in the limit i> 1 we get 




(3.39) 




(3.40) 



It can be real only if h = — tts (note that A( — 1) < 0), cf. (Macek et al. 2005). 
Finally, the result that we are interested in is the linear combination 



fe(x) = f(x) +jxe(x), 



(3.41) 



where the complex number 7 is chosen such that fe{x) contains only an outgoing wave 
at large x (this corresponds to an atom and a dimer flying apart after the three-body 
recombination event). This condition gives 



_. 7r^ A(0) v /27rA(-l) 



sinh(7rso) cosh(7rso — i6) 



Keeping only the relevant oscillating term at large x we obtain 



h{x) 



> i47sin0sinh(7rso) cxp(ia; + io~)/x. 



(3.43) 



X 



00 



So, we have found the coefficient in front of the outgoing atom-dimer wave, which 
is enough to calculate the atom-dimer outgoing flux. Indeed, the large- a; asymptote 
/ = £ exp(ix)/x, where £ is any complex amplitude, is accompanied by the flux |£| 5>ooi 
where 



$oo = 3 x (871-)- 1 x (4?r) x 2 = 3. (3.44) 

Here we have explicitly written out the following factors: the factor of 3 reflects the 
three symmetric possibilities of forming the dimer (corresponds to the interchange 
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ri r2 rs), the factor of (87r) _1 arises from the relation in between \I> and / 
(see the discussion after Eq. (|3.28p h the factor of 4tt is the solid angle in the outgoing 
atom-dimer channel, and the last factor of 2 is the atom-dimer relative velocity in the 
x, y-coordinates. The three-body recombination rate constant a s is obtained by taking 
the squared modulus of the prefactor in front of cxp(ix + io~)/x in Eq. ()3.43j) and by 
multiplying it by , by the factor of 1 /6 reflecting the fact that the number of triples 
in the gas is n 3 /3!, and by the factor fia 4 /m in order to restore the original physical 
units. We should also mention that the 9-dimensional volume V 3 is taken to be a 
unit volume in the original system of coordinates {ri,r2,r3}. In the new coordinates 
{x, y,R cm }, where R cm is the center-of-mass coordinate, this volume equals V 3 = 
8/3-\/3- The final result for the three-body recombination rate constant reads 



a s = 128^(4^-3^3) 



sin 2 [so ln(a/a*o)] + sinh 2 77* 



ha 4 



sinh (ttsq + 77*) + cos 2 [s ln(a/a *)] fn 



(3.45) 



where we have expressed the three-body parameter 9 through the original physical 
units: 

9 = s 111(0/0*0) + ir)*. (3.46) 

In Fig. 13.51 we plot the quantity mas/ha 4 as a function of a for different values of the 
elasticity parameter 77* . 




100 



Fig. 3.5 The rate constant a s for the recombination into a weakly bound state for three 
identical bosons. The function is log-periodic with the scaling factor exp(-7r/so) ~ 22.7. 



The rate of recombination into a weakly bound level is only a fraction of the total 
loss rate. The rest is due to the formation of deeply bound molecular states. We have 
already discussed this mechanism in Sec. 12.61 In contrast to the recombination into 
shallow states we now have to look at the balance of the incoming and outgoing fluxes 
of atoms corresponding to the short distance asymptote of fg (x) given by Eq. (|3.4ip 

f e {x) ► 7xe(x) = -f(Ae te x- 1+tS0 + A* e~ t6 x' 1 - 130 ), (3.47) 
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In analogy with let $0 denote the number of atom triples disappearing at the 
origin (x = 0, y = 0), provided the function / takes the form of the incoming wave 
x -i-is w jth unit weight. With this definition the recombination rate constant follows 
from Eq. (|3T47j) : 

a d = (l/3!)(fta 4 /m)| 7 | 2 $ |A| 2 2sinh(2^), (3.48) 

where the prefactor 1/3! is related to the number of triples in the gas. The product 
$o|^4| 2 can easily be found from the definition p.32[) by equating the fluxes at x — > 
and at x 00 and using Eq. (|3.44p : 

$ |^| 2 = 2$oo sinh(27rs ) = 6 sinh(27rs ). (3.49) 
Substituting Eqs. (|Q9|) and (f3T42|) into Eq. (pT48|) we obtain 

a d (a > 0) = 128^(4. - 3^3) coth^p) coshfa) sinhfa ) ^ 

sinh (-ksq + 77*) + cos 2 [so ln(a/a *)] m 

Since in the case of identical bosons the product 7rso is rather large the function 
ad{a > 0)/a 4 is almost flat (we can neglect the a-dependent term in the denominator) 
in contrast to a s /a A . The constant ctd(a > 0) monotonically increases with 77* and in 
the extreme limit 77* — > 00 the ratio between the two rate constants equals ad(a > 
0)/a s = coth(7Ts ) ~ 1.0036. 

Let us now discuss the negative side of the resonance and derive ay (a < 0). In this 
case we use the units H = \a\ = m = 1 and some equations described above should be 
modified accordingly. In particular, Eq. (|3.26[) reads 

(L + l)f(x) = 4tt*o (3.51) 

and we now introduce an auxiliary function g$ related to / by 

f(x) = 4tt*o[1 - A(0)/x + A(0)A(-l)ff O (x)], (3.52) 

where go satisfies (L + l)go(x) = x~ 2 . We write the solution in the form 

M*)=^rx(z)d Z r^ (3.53) 

27ra Jo Jo y + x 

and integrating it in the small- a; limit we get the asymptote 

(3 , 4) 

x->o smn(7rso) 

The function / is a solution of Eq. (|3.51[) . but its oscillations at small x do not have 
(in general) the correct phase imposed by Eq. (|3.33|) . This difficulty is resolved by 
observing that /* also satisfies Eq. (|3.51[) . The correctly behaving solution reads 

f 1 \ exp(ig)/(a;) +exp(-ifl)/*(a:) 

MX) = eM^+eM-m ' (3 ' 55) 
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and by subtracting the outgoing flux from the incoming one at small x we obtain the 
result 

coth(7rso) cosh(?7*) sinh(ry*) ha A 



a d (a < 0) = 1287r 2 (47r - 3^3)- 



COS' 



[soln(|o|/ao*)] + sinh (17*) m 



(3.56) 



For small 77, Eq. (|3.56|) is characterized by resonances at a = — ao* exp(-7r/2so) modulo 
exp(7r/so) (see Fig. 13. 6[) . These points mark the passage of Efimov trimers across the 
three-atom threshold with the usual consequences: the contribution of the trimer state 
(of size ~ \a\) in the three-body wavefunction becomes very large (close to 8 = ir/2 
the right hand side of Eq. ()3,55[) diverges) and atoms spend a lot of time close to each 
other, which leads to enhanced recombination losses. Note that on average a.d{a < 0) is 
significantly larger (by about three orders of magnitude) than a s and ay (a > 0) for the 
same \a\. In the limit 77* — > 00 we have otd(a < 0)/a s = cxp(27rso) coth(7rso) ~ 558.9. 
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Fig. 3.6 The rate constant ad for the recombination into deeply bound states for three 
identical bosons on the negative side of the Feshbach resonance. The peaks correspond to 
Efimov states crossing the three-atom threshold. 



Recalling that a_ is defined as the value of a where ay (a < 0)/a 4 reaches its 
maximum, Eq. (|3.56[) implies 

|a_|/a * = exp(7r/2s )- (3.57) 

The zero-range theory thus predicts that the maxima of a s and ay (a > 0) and the 
maxima of ad(a < 0) are placed at a = ±|o_|, i.e. symmetric with respect to the 
center of the Feshbach resonance. 

This analytical approach can be extended to homonuclear and heteronuclear mix- 
tures of fermions and bosons by taking into account the mass imbalance and/or nonzero 
angular momentum. It can be used as a zero-energy reference point for numerical 
approaches which are supposed to give results also at finite energies (three-body re- 
combination at finite temperatures, atom-dimcr scattering, etc.) It is also of immense 
help in the cases of small sq where the exponentially large scaling parameter impedes 
numerical calculations. 
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3.4 STM equation near a narrow resonance 

Let us now discuss how one can include the finite width of the resonance into the STM 
approach. In fact, the derivation of Sec. 13.11 is also valid in the narrow resonance case 
except that the STM equation itself ()3.14j) should be modified. Indeed, the asymptote 
(|3.12l) should be compared with the narrow resonance Bethe-Peierls condition (|1.17|) 
which in current notations reads 

tf(xi,yi) cx l/j/i - - y^lRlEx, y x -> 0, (3.58) 

where E\ is the collision energy of atoms 2 and 3, i.e. the kinetic energy of motion 
"along" the yi-direction of the six-dimensional {xi, yi}-space. Using Eq. p.8[) and 
the definition (|3. 15|) we write 

-V 2 * V 2 * V 2 fi 

El = lim ^- =E+ lim = E+ (3.59) 

yi^O yi ^o v]/ f 1 

Substituting Eq. ([339]) into Eq. ([338]) and comparing the latter with (|3.12|) we obtain 
the narrow resonance STM equation 

L ljB {/}( Xl )+ [V=E- - ^WiR\{E + V^)]/i(xi) = 47t* (xi, 0). (3.60) 

The equations for fa and f% should also be modified by including similar terms pro- 
portional to i?2 an( i ^3 j respectively. Equation (|3.60p can be, of course, written in 
momentum space. The operator — is then substituted by p\ and the correspond- 
ing representation of L\^e is given by Eq. (|3.16[) . 

At sufficiently short distances the operator V 2 X in Eq. (|3.60[) dominates and the 
short-distance asymptote of fi satisfies — V 2 1 /i(xi) = 0. This is consistent with 
Sec. 12.81 where we argued that at these distances the three-body wavcfunction de- 
scribes the free motion of an atom and a closed channel molecule. In fact, according 
to Eq. (|1.27[) the probability for atoms 2 and 3 to be in the closed channel equals 
4ttRI\ lim yi ^ 2/i * 1 2 = 47ri? 1 |/i(x :L )/47r| 2 , and, therefore, ^ R^/Anf 1 (x 1 ) can be con- 
sidered as the wavefunction of the relative motion of an atom and a bare molecule. In 
the narrow resonance limit they interact only by virtually breaking up the molecule 
and exchanging one of its constituents with the free atom. This exchange interaction 
is given by the integral operator L\_e in Eq. (|3.60j) and can be treated pcrturbatively 
in the limit R* — > oo. Obviously, we can also introduce the direct interaction between 
the closed channel molecule and atom by imposing another zero-range Bethe-Peierls 
boundary condition, this time on the function f\. 

3.5 Atom-dimer scattering near a narrow resonance 

To illustrate how the approach of Sec. 13.41 can be implemented in practice let us 
consider a concrete example. Namely, we calculate the atom-dimer scattering length 
for a system of two statistically identical atoms of mass M and another atom of 
mass m. We use the notations of Sec. 13.21 with upper sign for identical bosons and 
lower - for fermions. To account for the finite resonance width we now add the term 
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-y/2jlR*(E + V 2 )/(x) to the left hand side of Eq. (|3T7]) and look for its solution at 
zero atom-molecule collision energy, i.e. we set E = eo < 0, where eo is the energy of 
the molecular state [see Eq. (|1.31[) ]. Since E is negative, the free solution Vfo vanishes, 
and Eq. (|3.17|1 takes the form 

{-^iR*Vi + L B=eo )/(x) =0. (3.61) 

The algorithm of calculating the atom-dimer scattering length a a d is straightfor- 
ward: expand Eq. (|3.61l) in spherical harmonics, solve it for I = (s-wave symmetry), 
and deduce a a d from the large- a; asymptote of f(x). Indeed, at distances x ^> 1/k, 
where k = y/2fx\eo\, Eq. (|3.61[) describes the free atom-molecule motion, i.e. the s-wave 
symmetric solution is a linear superposition of x~ 1 and x . The relation between the 
corresponding coefheicnts is fixed by the atom-dimer scattering length: 



f(x) > N(l - y / 2pa ad /x), (3.62) 

x— >oo 

where p = m(M + m)/(2M + m) is the atom-molecule reduced mass and N is the 
normalization prcfactor, which is actually not needed for determining a ac i, but we will 
return to it when discussing the inelastic atom-dimer relaxation. 

Alternatively, one can follow the same procedure in momentum space. Introducing 
the Fourier transform /(p) = / /(x) exp(— ipx)d 3 x we look for the solution in the 
form 

f(p) = N[(2tt) 3 S( P ) + 4irg(p/V^)/(^>p 2 )}- (3.63) 

We substitute p.63[) into the momentum space STM equation, integrate over the 
angles of p and arrive at the following equation for the function <?(fc), 

^2^R*g{k) = ±- 



(fc 2 + COS 2 (f) COS ( 



VT+W-l .... 1 rg(k'). fc 2 + fc' 2 + 2fcfc'sin^ + cos 2 <^ , 

-9( k ) ± ■ / TTT lo § 75 , 1,9 I I —3 A k > ( 3 - 64 ) 



k 2 * w 7T sm(2(/)) J kk 1 b k 2 + k' 2 - 2kk' sin <j> + cos 2 , 



the solution of which gives the atom-dimer scattering length: a ac i = — g(0)/\/— 2fleQ. 

In Fig. 13.71 we plot a a d in units of a calculated from Eq. (|3.64p for collisions of 
K atoms with K-Li molecules formed near a narrow interspecies resonance. The solid 
line shows the case of identical bosonic 39 K atoms and the dashed one - identical 
fermionic 40 K. We see that in the limit R*/a <?C 1 the bosonic case is Efimovian and 
is characterized by the log-periodic dependence of aw/a on a of the form (|2.22p . the 
parameter R* plays the role of the three-body parameter as we have discussed in 
Sec. 12.81 In contrast, in the fermionic case the atom-dimer scattering length has a 
well-defined limit for R* — and the left hand side of Eq. (|3.64[) can be considered as 
a weak perturbation for small R* /a. 

In the opposite limit, R* /a^$> 1, which we call the regime of intermediate detuning, 
we treat Eq. (|3.64p pcrturbatively in the small parameter r\ = 2/ieo-R* <C 1. 

Namely, we look for the solution in the form g(k) = g(°'(k) + g^ x \k) + g^ 2 \k) + 
where the functions g^ are of order rf and can be found recursively: g^ = 0, g^ 1 ' (k) = 
±7?cos _1 (0)/(fc 2 +cos 2 0), then <?W for i > 1 is obtained by substituting g' 1-1 ' into 
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Fig. 3.7 The atom-dimer s-wave scattering length for 39 K-( 39 K- 6 Li) (solid line) and 
40 K-( 40 K- 6 Li) (dashed line) collisions near a narrow interspecies resonance. The dotted lines 
show the large R* /a asymptotes (|3.67|l . 



the <?-dependent terms in the right hand side of Eq. (|3.64[) and by multiplying the 
result by r\. In this way we obtain the first two leading terms 



?(0)«±-5f 



— 1 ± cos 2 (b 

1 + o -n 



(3.65) 



Then, expressing the small parameter in terms of \J a/R* and using 



' = vi + ^/.-i ~V^( 1 + 5V^) <3 - 66) 

wc finally obtain the atom-dimer scattering length up to the next-to-leading order in 

— * T oo^ \ l + 2 V ^) ■ ( } 

The dotted lines in Fig. 13.71 show the asymptotes p.67|) . 

Let us now discuss how the atom-dimer relaxation can be calculated for molecules 
formed near a narrow resonance. The process is local, it happens at distances of the or- 
der of the size of the closed channel molecule, which is much smaller than R* according 
to our definition of narrow resonances. We have already mentioned that the probability 
of finding three-atoms in a small volume of size -C R* is dominated by the probabil- 
ity of finding there an atom and a closed channel molecule. Therefore, the probability 
density for finding three atoms in the recombination region equals (R* / 4ir)\ f (x = 0)| 2 , 
which is proportional to the relaxation rate if we treat the process perturbatively. The 
proportionality prefactor requires solving the three-body problem at short distances, 
which is a challenging task. However, since the shape of the closed channel wave- 
function is not dramatically sensitive to the magnetic field, we can assume that this 
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microscopic prcfactor stays approximately constant close to a given Feshbach reso- 
nance. Therefore, the quantity that can be calculated in the zero-range approximation 
is the ratio a re i(a)/cvbare of the relaxation rate constant at a given a to the relaxation 
rate constant for collisions of atoms and closed channel molecules. This ratio should 
tend to 1 in the limit R* ^s> a. 

In Fig. 13.81 we plot the ratio a re i(a)/abare for the same physical systems and us- 
ing the same notations as in Fig. 13.71 i.e. the solid line stands for bosonic 39 K and 
dashed - for fermionic 40 K. The results are obtained by substituting the already cal- 
culated function g(k) into Eq. (|3.63[) and evaluating f(x = 0) by integrating f(k) over 
momenta. One should keep track of the a-dependence of the normalization prcfactor 
N oc anorm given by Eq. ()1.35j) . This dependence follows from the correct normaliza- 
tion of the dimer wavefunction in the narrow resonance case (see Sec. I1.3[) . The final 
result reads 

2 



a re i(a) _ y/1 + 4R*/a - 1 
a 



1 + - 



g(k)dk 



(3.68) 



The dotted lines in Fig. 13.81 show the large R* /a-asymptotes and are obtained by 
approximating g(k) w g^\k) in Eq. (|3.68[) . 
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Fig. 3.8 The atom-dimer s-wave relaxation rate constant versus the detuning. The physical 
systems and notations are the same as in Fig. 13.71 The dash-dotted lines reflect the power 
scaling of a rc i at large a (see text). 



In the regime of small detuning, R* /a <C 1, the dependence of cvrei/^bare on R* ja 
can be estimated qualitatively. Consider first the fermionic case. The function f(x) is 
of order N ~ 1/y/a at distances ~ a where it should be matched with the power law 
f(x) ~ a~ 1 / 2 (x/a,y i approximately valid in the window R* < x < a. The power v\ is 
the root of \i =n {v) [sec Eq. (j3~22|) and Fig. GO]. For the system 40 K- 40 K- 6 Li it equals 
v\ w 2.0193. Then at distances x < R* the function / is approximately constant since 
the interaction in between the atom and the closed channel molecule is not resonant. 
We finally get 
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a r ci/abarc oc R*\f(x = 0)| 2 oc (R* /a)- 2ui -\ R* < a. (3.69) 

This qualitative scaling is shown in Fig. 13.81 as the lower dash-dotted line. 

We see that the s-wave atom-dimer relaxation is highly suppressed in the limit 
R* <C a. A good news is that even for a = R* , where the qualitative scaling (|3.69p 
shows no suppression, the quantitative result predicts the suppression of the relaxation 
by three orders of magnitude compared to the "bare" atom-molecule case. Apparently, 
for this relatively large mass ratio the exchange of the light atom leads to an appreciable 
repulsion between the heavy ones (remember that according to the Born-Oppenhcimer 
approach, the light atom is in the "repulsive" ip- state). These effects are not taken 
into account in the derivation of Eq. (|3.69p . 

As far as the bosonic case is concerned we can repeat the same speculations disre- 
garding the atom-dimer resonances and the oscillatory part of f(x) ~ a~ 1 ^ 2 x~ 1 cos[so ln(x/ro)] 
at distances R* < x < a. We obtain then for bosons a re i/abarc ~ (-R*/a) -1 ! i-e. 
Eq. (|3.69[) with V\ = — 1 (the upper dash-dotted line in Fig. I3.8|) . Figure l3~8l clearly 
shows that in this case there is no suppression of the relaxation. 
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Final remarks 



4.1 Outlook 

Let us now briefly mention other rapidly developing directions of the theoretical few- 
atom physics. Molecular regimes in fermionic and bosonic mixtures require a better 
understanding of collisional properties of molecules, which is a four-atom problem. The 
transition from three-body to four-body is conceptually straightforward, it requires 
just a bit more space for formulas, although much more computing power. The four- 
body STM equation for fermions in coordinate and momentum space has been derived 
along the lines of Chap. H by Petrov et al. (2004, 2005a, 2005b). The three-body 
STM approach has also been generalized to the four-body case by using diagrammatic 
techniques (Brodsky ct al. 2005; Levinscn and Guraric 2006). The Born-Oppcnhcimer 
approximation has been applied for calculating collisional properties of highly mass 
imbalanced heteronuclear molecules (Marcclis et al. 2008). Various aspects of the four- 
boson problem with short-range interactions (universality, tctramer states, necessity 
of an additional four-body parameter) have been actively studied in momentum space 
by using the Faddeev-Yakubovsky equations (Platter et al. 2004a; Yamashita et al. 
2006, Hammer and Platter 2007; Deltuva 2010), by using the adiabatic hyperspherical 
approach (von Stecher et al. 2009; see also review of Rittenhouse et al. 2011), and by 
other methods (Hanna and Blume 2006; von Stecher 2010; Yamashita et al. 2010). The 
problem of three heavy fermions interacting with a light one has recently attracted 
attention (Castin et al. 2010; Blume and Daily 2010; Gandolfi and Carlson 2010; Mora 
et al. 20115). In particular, Castin et al. (2010) have argued that in this system the 
four-body Ehmov effect can occur for mass ratios M/m < 13.6, i.e. in the absence of 
the three-body one. 

Another direction of research, obviously relevant for ultracold gases, is the few- 
body problem in a trap. The spectrum of three trapped bosons has been calculated 
by Jonsell et al. (2002). The relation of the trapped problem at unitarity (a = oo) 
and the free-space problem at zero energy has been established by Tan (2004) and by 
Werner and Castin (2006). Various results on the spectrum and correlation functions 
of few-fermion systems in a harmonic potential have been obtained by Blume and 
co-workers (von Stecher et al. 2007, 2008; Blume et al. 2007). For an overview of the 
few-body problem in a trap see (Blume 2012). 

A notable progress has been made in solving the few-body problem in low dimen- 
sions. The hyperspherical approach has been applied to the two-dimensional three- 
boson problem by Nielsen ct al. (1999) and by Kartavtsev and Malykh (2006). Platter 
et al. (2004b) have calculated the bound state energies of four two-dimensional bosons 
by using the Faddeev-Yakubovsky equations. Brodsky ct al. (2005) have considered 



42 Final remarks 



the two-dimensional few-body problem in the diagrammatic framework. Heteronu- 
clear two-dimensional trimers have been discussed by Pricoupcnko and Pedri (2010) 
by using the two-dimensional STM equation in momentum space. Liu et al. (2010b) 
have calculated the second and third virial coefficients for a two-dimensional strongly 
interacting Fermi gas. 

As far as the one-dimensional case is concerned the problem of N particles is exactly 
solvable in some cases (Lieb and Liniger 1963; Lieb 1963; McGuire 1965; Yang 1967; 
Gaudin 1967) and the corresponding wavefunctions can be written out explicitly. The 
exact results can be compared to the results of various few-body methods (Dodd 
1970; Thacker 1975; Amaya-Tapia et al. 2004). Speaking of nonintegrable cases, the 
atom-dimer and dimer-dimer scattering problem in quasi-one-dimensional geometry 
has been studied by Mora and co-workers (Mora et al. 2004. 2005a. 2005b). Muryshev 
et al. (2002), Sinha et al. (2006), and Mazets et al. (2008) have perturbatively estimated 
the three-body diffraction effect originating from the finite transversal size in a quasi- 
1D system of bosons. Deviations from integrability originating from the finite width of 
a Feshbach resonance have been investigated by Yurovsky et al. (2006). The spectral 
statistics and the response of ID few-body systems with and without mass imbalance 
(nonintegrable and integrable, respectively) have been studied by Colome-Tatche and 
Petrov (2011). The one-dimensional three-body problem on a lattice has recently been 
discussed by Orso et al. (2011) and Valiente et al. (2010), see also (Keilmann et al. 
2009). 

Different components of a heteronuclear mixture can feel different external forces. 
For example, we can imagine that one of the components is confined to a quasi- 
one-dimensional or quasi-two-dimensional geometry whereas the other component re- 
mains free. Besides, by adding an optical lattice one can separately control the effec- 
tive masses of the components. This leads to interesting peculiarities of many-, few-, 
and even two-body problems in such mixed-dimensional systems (Petrov et al. 2007; 
Nishida and Tan 2008; Levinsen et al. 2009; Lamporesi et al. 2010). 

We should mention the analysis of the few-fermion problem near a p-wave resonance 
(Levinsen et al. 2007; Jona-Lasinio et al. 2008) which is interesting in view of the 
realization of nontrivial superfluid phases in strongly interacting polarized fermionic 
gases. Finally, the progress in preparing cold dipolar gases (see the lecture course of 
D. S. Jin and J. Ye in this book) motivates studies of the few-body problem with 
dipolar interactions (Wang et al. 2011b, 2011c). 
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